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The  study  concludes  that  yes,  there  are  physically  realistic  para- 
meters that  can  be  found  (using  Parameter  Estimation  Techniques)  that 
allow  slabs  to  be  adequately  modeled  this  way  and  that  reproduce  that 
actual  dynamic  response  of  slab.  Secondly,  these  parameters  can  often 
(but  apparently  not  always)  be  determined  from  static  test  data,  elemen- 
tary theory,  or  combinations  thereof . Recommendations  are  made  concerning 


the  future  design  and  testing  of  R/C  slabs. 
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1. 


INTRODUCTION 


Estimating  the  vulnerability  of  a structure  to  an  attack  by  a nuclear 
weapon  requires  some  method  of  determining  the  response  of  the  structure 
to  the  transient  loadings  produced  by  the  detonation  of  the  weapon.  Single 
degree-of-freedom  systems  are  commonly  used  within  the  defense  community 
as  the  means  of  modeling  the  response  of  structures  to  these  transient 
loadings.  While  it  is  recognized  that  estimates  of  strength  and  stiff- 
ness produced  by  single-degree-of-f reedom  system  models  may  be  in  error, 
it  has  been  generally  argued  that  the  models  are  adequate  since  the  over- 
alluncertainties of  the  problem  do  not  justify  the  use  of  more  complex 
system  models.  The  objective  of  this  study  is  to  exmamine  the  question  of 
the  adequacy  of  a single-degree-of-freedom  model  to  predict  the  transient 
response  of  one  class  of  structural  member,  the  slab.  The  approach  used 
is  to  seek  the  answer  to  two  related  questions: 

1.  Is  there  any  set  of  parameters  for  a single  degree-of-freedom 
model  that  will  match  the  experimentally  derived  transient  re- 
sponse of  various  slab  designs? 

2.  Can  these  parameters  reasonable  be  estimated  on  an  a priori 
basis  from  either  static  test  data  or  from  engineering  principles 

The  analytical  procedure  used  is  the  technique  of  Parameter  Estima- 
tion. This  technique  basically  takes  a system  model  of  any  degree  of 
complexity  and  determines  that  set  of  model  parameters  which  minimises 
the  differences  between  the  observed  and  calculated  response  values. 

The  report  itself  is  organized  in  the  following  manner.  The  ideal- 
ization of  the  slab  as  a single  degree-of-freedom  system  is  discussed 
rirst,  (i.e.,  the  assumptions  and  approximations  made),  followed  by  a 
discussion  of  the  "multi-parameter"  one  degree-of-freedom  model  used 
in  the  study.  The  next  section  deals  with  estimation  of  parameters  (i.e., 
fitting  the  model  to  the  data) , followed  by  results  presented  for  three 
(widely  different  R/C  slabs.  A section  of  suggested  procedures  for  ana- 
lyzing slabs  as  single  degree-of-freedom  models  is  then  presented,  and 
the  report  ends  with  conclusions  and  recommendations  for  the  future  design 
and  testing  of  such  structures. 


IDEALIZATION  OF  A SLAB  AS  A SINGLE  DEGREE-OF-FREEDOM  SYSTEM 


2. 

This  section  follows  the  discussion  by  Watt  (Reference  3) , and  simi- 
lar results  are  also  presented  in  Biggs  (Reference  2)  and  other  texts 
(Reference  4).  In  the  design  of  blast-resistant  structures,  considerable 
economy  can  be  realized  if  the  design  includes  not  only  the  elastic  but 
also  the  plastic  resistance  of  the  structure.  Plastic  behavior  is  not 
generally  permissible  under  continuous  operating  conditions,  but  is  in 
many  cases  quite  appropriate  for  design  when  the  structure  is  subject  to 
a severe  dynamic  loading  only  once  or  twice  during  its  life. 

It  can  be  concluded  from  examination  of  exact  or  rigorous  dynamic 
analysis  (i.g.,  Reference  2)  that  relatively  simple  structures  having 
simple  boundary  conditions  can  be  easily  analyzed.  Otherwise, the  analysis, 
though  not  impossible,  becomes  cumbersome.  For  this  reason  and  for  practi- 
cal design  purposes,  the  idealized  spring-mass,  single-degree-of-f reedom 
approximate  analysis  is  commonly  used.  This  method  provides  a means  for 
analyzing  the  more  complex  structures  rapidly  and  with  a reasonable  de- 
gree of  accuracy  (Reference  2). 


2.1  EQUATION  OF  MOTION 


The  actual  deep  slab  structure  and  the  idealized  spring-mass  system 
are  shown  in  Figure  2-1.  The  first  peak  value  of  deflection  is  of  primary 
interest  for  this  problem,  and  damping  is  often  neglected  (Reference  4), 
but  for  completeness  it  has  been  retained.  From  the  free-body  diagram  shown 
in  Figure  2-1,  the  equation  of  motion  for  the  equivalent  1 d-o-f  system  is 


where 


y 


y 


F 


e 


Mey(t)  + 6ey(t)  + Key(t)  “ Fe(t)  = 0 (2_1 

2 

Me  = equivalent  or  effective  mass,  lb-sec  /in  (Kg) 

(t)  = acceleration  of  slab  at  time  t,  in/sec2  (m/sec2) 

Be  = equivalent  or  effective  damping,  lb-sec/in  (Newtons-meters/ 

sec) 

(t)  = velocity  of  slab  at  time  t,  in/sec  (m/s ec) 

Ke  = equivalent  or  effective  stiffness,  lb/in  (Newtons/m) 

(t)  = deflection  of  slab  at  time  t,  inches  (m) 

(t)  = equivalent  or  effective  force  at  *' ■■■  ■■-nnuti  . 

2-1  ; PR  f?CEoiY\V  ^ C" 


AIRBi-AST  pressure,  P(  t ) 
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Figure  2-1.  Deep  Slab  Specimen  and  1 d-o-f  Model 


In  order  tc  define  the  equivalent  system,  the  parameters  Mg,  g^ , K , 
and  Fe  must  be  evaluated.  The  equivalent  one-degree-of  freedom  system  is 
one  for  which  the  kinetic  energy,  dissipation,  internal  energy,  and  work 
done  by  all  external  forces  are  at  all  times  equal  to  the  same  quantities 
for  the  continuous-mass  deep  slab  structure.  It  is  assumed  that  the  de- 
flected shape  of  the  structure  is  the  same  as  that  due  to  static  loading. 
(This  is  a key  assumption.) 

If  any  point  on  the  deflecting  surface  of  the  deep  slab  is  described 
by 


y (x,  z , t)  = A (t)  0 (x,  z) 


(2-2) 


where 


A ( t ) = displacement,  inches  (meters)  as  a function  of  time  t 
0(x,z)  = assumed  deflecting  shape  (i.e.,  mode  shape)  of  the  deep  slab 

the  velocity  y(t)  of  any  point  on  the  slab  becomes 

y (x , z , t)  = A(t)0(x,z)  (2-3) 

If  the  equivalent  system  is  to  respond  similarly  to  the  midspan  of 
the  actual  structure,  then  the  displacements  and  velocities  of  both  systems 
must  be  the  same. 


Equating  the  kinetic  energy  of  both  systems  yields 
l/2MeA(t)2  = 1/2  fS  M [ A ( t ) 0 (x,z)  ] 2 dS 


or 


-/ 


M0  (x,s)  dS 


(2-4) 


where  S is  the  area  of  the  slab  (planform) . 


Equating  external  work  yields 


FA  (t) 


/ 


F (A  (t)  0 (x,  z)  ] dS 


-r 


F 0 (x, z)  dS 


(2-5) 


2-3 


L 


where 


2,3  2 

M = mass  per  unit  area,  lb-sec  /in  (kg/m  ) 

2 2 

F = force  per  unit  area,  pounds/in  (Newtons/m  ) 

Similarly,  equating  the  dissipation  yields 
rS  . 

BeA(t)  =J  cA(t)*(x,z)  dS 

where  c is  a "localized"  damping  term. 

The  equivalent  load  of  Equation  (2-5)  applies  to  the  magnitude  of  the 
load,  while  both  loads  have  the  same  time  function. 


The  resistance  or  stiffness  is  the  internal  force  tending  to  restore 
the  system  to  its  unloaded  static  position.  If  the  stiffness  is  defined  in 
terms  of  the  load  distribution,  then  the  maximum  stiffness  is  numerically 
equal  to  the  total  load  of  the  same  distribution  tp(x,z)  which  would  cause  a 
unit  deflection  at  the  point  where  the  deflection  is  equal  to  that  of  the 
equivalent  system. 


The  shape  function  <p(x,z)  changes  as  the  deep  slab  progresses  through 
the  different  stress  ranges,  i.e.,  elastic,  elastic-plastic,  and  plastic. 
Thus,  a complete  solution  requires  that  a shape  function  for  each  range  be 
considered.  Or,  to  put  it  another  way,  the  effective  mass,  Mg,  is  really 
Me(y),  i.e.,  it  is  displacement-dependent . Similarly,  6e  and  Fg  depend  upon 
the  displacement. 

2.2  CANONICAL  FORM  OF  THE  EQUATION  FOR  VARIOUS  MODE  SHAPES 


If  one  defines  the  total  mass  (of  the  slab)  by 

Mt  =J  Ph  dS 


then 


Mt  = phS 


( 2 — 6 a ) 


( 2— 6b ) 


where  S is  the  planform  area  of  the  slab  and  it  is  assumed  that  the  average 
density  (p)  and  the  slab  thickness  (h)  are  constant.  Similarly,  the  total 
force  applied  is  defined  by 


and  when  the  pressure,  p(t)  is  assumed  to  be  uniform  (over  the  planform 
area)  one  has 


1 ( 2- 7b) 

Finally,  if  one  combines  the  preceding  results  (Equations  2-6b,  2-7b, 

and  2-4,  2-5,  with  2-1)  and  then  divides  through  by  the  area,  S,  he  obtains ’ 
the  canonical  form 


- MY  + 6gy  + ke(y)  = ap(t) 
where  the  mass  coefficient 

u _ li^Ph  Q2  dS 
[ chs 

and  the  u is  always  1 1.  Similarly  kg  (y)  is  the  restoring  force 

(per  unit  area)  (i.e.,  it  represents  the  static  load-deflection  curve, 
pressure  vs.  deflection).  it  is  clear  that 


(2-8) 


(2-9a) 


J/xi 

( PS 


(J5  p(t)  *■ 


(2-9b) 


and  that  a < 1. 

Estimates  of  the  coefficients  U and  a can  be  obtained  by  assuming 
various  forms  for  the  mode  shape,  0.  For  example,  when  0=1,  (i.e.,  the 
slab  translates  through  space  like  a rigid-body),  then  p=l  and  a=l. 

xlarly,  if  . - sin  — (bending  like  a simply-supported  beam  in  one 
direction  only),  then 


. ./si"2v 


/•*»  T 


' (I) 


Or,  if  the  slab  acts  like  a plate  (or  membrane)  in  two  directions,  then 


<Hx,z)  = sin  — sin  ~ 
a b 

u = i bffsin2(ir ) siT[2($  dxdz  = (§)  = -25 

“ = ib/J (sin  ir  sin  tJ1) dxdz  = (l)  = -4053 

Thus,  a priori,  one  has  bounds  on  the  coefficients  a and  p of  the 
canonical  Equation  (2-8).  These  "bounding  values"  will  be  used  subsequently 
when  the  parameters  of  Equation  (2-8)  are  estimated.  As  stated  previously, 
kg(y)  represents  the  static  load-deflection  curve,  and  B y is  the  damping 
term,  which  has  practically  negligible  influence  on  the  maximum  initial 
response  (cf . Biggs,  Reference  2) . A typical  pressure  vs.  deflection  curve 
is  shown  in  Figure  2-2,  and  Figure  2-3  shows  a representative  deflection- 
time history. 

2.3  "REAL-WORLD"  EXPERIMENTAL  EFFECTS 

In  practice,  the  mass  coefficient  p (and  the  force  coefficient,  a) 
will  vary  with  the  displacement  y.  For  example,  a thick  slab  might  be  ex- 
pected to  start  out  deflecting  linearly  (with  p = .5,  say),  and  end  up 
"punching  through"  (with  p = 1,  say)  as  the  pressure  loading  increases.  In 
other  words,  one  might  have 

. 5 £ p (y)  <1  (2-10a) 

during  the  dynamic  response.  Similar  restrictions  apply  to  the  force  coeffi- 
cient, a, 

. g < a (y)  < 1 ( 2-10b) 

for  thick  slabs.  This  type  of  behavior  was  observed  by  Watt,  Reference  3. 

Conversely,  for  thin,  conventional  slabs,  Brown  and  Black  (Reference 
5)  found  that  they  began  deflecting  linearly  (like  a plate)  and  ended  up 
stretching,  like  a membrane.  In  this  case,  one  might  expect  (at  least 
tentatively) 

.25  < p < .35  (2-] la) 

2-6 
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000* 


pical  Deflection  Tine-History  (Keenan  Reference  6) 


and 


. 4 I a 1 . 5 


where  p = p (y)  and  a = a (y)  both  vary  somewhat  with  deflection. 


Finally,  Keenan's  laced  slabs  (Reference  6)  behave  much  like  a wide 
beam  with  plastic  hinges.  Again,  one  might  estimate 


.25  < p < .33 


(2-12a) 


.51  a 1.6 


( 2- 1 2b) 


2 . 4 TEXTBOOK  RESULTS 

Biggs  (Reference  2)  presents  approximate  methods  for  the  design  of 

beams  and  slabs.  In  Chapter  Five  (pp.  199-244),  Biggs  discusses  a "mass 

factor"  Km  (analogous  to  the  p used  herein)  and  a "load  factor"  ^(analogous 

to  a).  Tabulated  values  of  K (p)  and  KT  (a) are  given  therein  for  beams  and 

m L 

slabs  subjected  to  various  dynamic  loads,  including  uniform  pressure. 

An  important  observation  made  by  Biggs  (and  not  originally  recognized 

by  the  authors)  is  that  the  factor  a (or  KT ) is  inherently  contained  in  the 



static-load  deflection  curve.  Thus,  the  equation  of  motion  can  be  written 


pMy  + By  +ak(y)  = ap( t) 


(2-13) 


or,  (dividing  through  by  a) 


MAy  + y + k (y)  = p(t) 


where  A=—  is  called  "the  load  mass  factor" 
a 


(2-14) 


In  equation  (2-14)  k (y)  is  the  static  resistance  function  (pressure 
vs.  deflection)  of  the  slab.  If  we  neglect  damping  (as  Biggs  suggests)  then 
equation  (2-14)  becomes 


MAy  + k (y)  = p(t) 


(2-15) 


Finally,  if  k(y)  is  approximated  by  a bi-1 inear  (elasto-plastic ) curve , 
k(y)  requires  just  two  parameters  to  define  it,  namely  piin(the  yield  stress) 
and  dlin  (the  displacement  when  yielding  begins).  Thus,  equation  (2-15)  has 


just  three  parameters  (X,  Plin,  and  dlin)  as  opposed  to  the  more  complicated 
multi -parameter  model,  which  will  be  discussed  shortly. 

Initially,  good  results  were  obtained  ( in  fitting  the  dynamic  response 
data)  using  the  multi-parameter  model.  However,  as  more  dynamic  data  were 
analyzed,  it  became  clear  that  the  most  physically  significant  results  would 
come  from  as  simple  a model  as  possible.  For  this  reason,  equation  (2-15) 
was  eventually  used  in  attempting  to  fit  the  dynamic  data. 


MULT I -PARAMETER  1 d-o-f  MODEL 


3 . 

When  this  study  was  begun,  it  was  not  known  how  well  a "simple  model" 
(with  constant  mass  p and  force  coefficient,  a)  would  perform  in  trying  to 
approximate  the  dynamic  response  data.  Hence,  the  one  d-o-f  model  was  made 
fairly  general  in  the  computer  program  (PEBLS) * which  was  used  in  this  work. 
If  one  rewrites  the  equation  of  motion  (2-8)  (repeated  below  for  convenience) 

uMy  + 6gy  + k (y)  = aP(t)  (2-8) 

and  notes  that  the  mass  term  p = p(y)  and  the  force  coefficient  a = a(y) , 
then  one  has 


Mp(y)y  + Bey  + kg(y)  = a(y)P(t) 


(3-1) 


Based  upon  the  limits  expected  for  (a,p)  (see  Section  2.  ) then  one 
might  choose  a linear  variation  of  (p,a)  with  deflection  as  shown  in  Figure 
3-1. 


Figure  3-la.  Variation  of  Force  Coefficient  » with  Deflection 


♦Parameter  Estimation  for  Blast-Loaded  Structures 
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If  the  pressure,  p(t),  is  measured  incorrectly,  (i.e.,  the  dynamic 
pressure  gage  is  not  reading  properly)  then  one  might  want  to  "scale  up" 
(or  scale  down)  the  pressure,  p.  Thus,  PEELS  uses 


P(t)  = Pscalef(t) 


I kj  1 


where  f (t)  is  the  experimental  (pressure  vs.  time)  data.  Finally,  the 
damping  term  w< 
little  effect. 


damping  term  was  kept  linear  (i.e.,  just  f J)  since  it  was  expected  to  have 


Thus,  for  generality,  the  eleven  parameters 
("o'  ^1'  ^d^  ~ rnass  coefficient 
(a.  a,,  a.)  ~ force  coefficient 

O j.  G 

(plin'  dlin ' Eplast)  ~ *i-linear  spring 

p . ~ pressure  scaling 

'scale  r ^ 

B ~ damping 

were  incorporated  in  the  computer  program.  The  resulting  model  is  referred 
to  herein  as  a "multi-parameter  1 d-o-f  model." 

Note,  of  course,  that  the  model  can  be  readily  simplified.  For 
example,  by  setting 

= C (or  ad  = 0) 

one  obtains 

, = = constant  (or  a = constant) 

i.e.,  one  can  revert  back  to  the  "constant  mode"  model  with  non-varving 
mass  anc  force  coefficients. 

Similarly,  by  putting  Epiast  = .001*  (or  .01)  one  can  achieve  an 
elastic  - perfectly  plastic  force-deflection  curve  for  the  restoring  spring. 


*Do  not  use  E_,  . = 0 or  B z 0 in  the  computer  program,  since  they  appear 

p 1 3 S XL 

in  the  denominator  within  the  F0RTRAN  coding. 

3-3 





. A*.  — ■ - Wt  ■ .■*. 


. 001* 


I 


Putting  Pscale  = 1-0  gives 
gives  low  damping.  Thus, 
multi-parameter  model  with 
just  noted. 


no  magnification  on  the  pressure,  and  £ = 
at  the  outset,  it  was  desirable  to  program  a 
the  flexibility  of  progressive  simplification 


Having  chosen  a "multi-parameter"  1 d-o-f  model,  one  is  then  faced 
with  the  problem  of  selecting  its  parameters:  ao,  , cx^...  etc.  In  most 
cases,  one  has  theoretical  results  to  guide  him,  (e.g.  C <_  u <_  1)  , and 
he  may  also  have  some  estimate  of  the  mode  shape,  d, , {which  allows  him 
to  compute  t , a)  and  in  many  cases  a static  load-def lection  curve.  With 
these  a priori  estimates,  the  analyst  has  a set  of  "initial  parameters" 
for  his  model.  Based  upon  his  experience,  the  analyst  also  places  a quanti- 
tative measure  of  confidence  on  his  initial  estimates,  and  these  numerical 
confidence  estimates  are  used  in  the  search  for  a final  set  of  parameters. 
This  search  for  "best  estimates"  of  the  parameters  is  termed  "parameter 
estimation,"  as  discussed  in  the  next  section. 


L 


*See  footnote  previous  page. 
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4. 


FITTING  THE  MODEL  TO  THE  DATA:  PARAMETER  ESTIMATION 


4.1  INTRODUCTION  AND  BACKGROUND 

The  question  of  selecting  the  parameters  (aQ,  etc.) 

of  the  1 d-o-f  model  in  such  a way  that  the  analytical  model  "matches"  the 
experimental  data  is  a problem  of  "parameter  estimation."  Parameter  estima- 
tion is  related  to  the  more  general  problem  of  "system  identification" 
(References  7 and  8).  The  main  idea  is  to  use  differences  between  measured 
and  predicted  behavior  (e.g.,  the  dynamic  response,  y,  as  a function  of 
time,  t)  to  adjust  key  parameters  of  a model  automatically  so  as  to  mini- 
mize response  differences  (i.e.,  minimize  y . , . , - y ,). 

-'analytical  -'experimental 

The  procedure  is  quite  general  and  can  be  applied  to  either  linear  or  non- 
linear mathematical  models  (cf.  Figure  4-1). 


INPUT 


i 1 


INITIAL 

PARAMETER 

VALUES 


Figure  4-1.  Flow  Diagram  of  Model  Verification/Parameter 
Estimation  Procedure 


In  the  present  study,  the  model  parameters  are  chosen  to  minimize 


the  vector  norm 


(4-1) 


N = { Au}T{Au} 
where  the  vector  {Au} 


U1  ^model 
Au2  = ^model 


is  given 

^expt^ 

■^expt' 


by 

at 

at 


time 

time 


t 

t 


1 

2 


“■ui  ^model 


-^expt^ 


at  time  t 


i' 


and  so  forth. 


Note  that  uncertainties  must  be  recognized  in  both  the  analysis 
(model)  and  the  test  (experiment).  That  is,  one  expects  (from  a static 
load-deflection  curve,  say)  to  know  the  yield  force  (stress)  fairly  well, 
but  there  is  still  an  uncertainty  in  it  (e.g.,  ±10%).  Similarly,  in 
the  measured  response,  Vexptf  one  does  not  know  y exactly , but  rather 
with  some  uncertainty.  Computational  methods  are  available  for  utilizing 
this  information  in  a quantitative  manner  so  as  to  achieve  an  optimum 
(minimum  variance)  fit  between  analysis  and  test.  In  other  words,  when 
confidence  in  the  data  is  greater  than  confidence  in  the  model  (analysis), 
the  model  is  adjusted  so  that  it  tends  to  match  the  data.  When  the  converse 
is  true,  changes  to  the  model  will  be  relatively  small. 


4.2 


STATISTICAL  PARAMETER  ESTIMATION 


As  just  discussed,  methods  for  parameter  estimation  which  recognize 
the  respective  uncertainties  in  both  analysis  and  test  have  an  obvious 
practical  advantage.  Such  a method  was  employed  in  this  study.  The  origi- 
nal derivation  and  application  of  the  method  appear  in  reference  9. 

The  particular  algorithm  for  statistical  parameter  estimation  des- 
cribed in  that  paper  uses  a linear  estimator  in  that  it  operates  on  the 
basis  of  a linearized  relationship  between  a parameter  vector,  (r),  and  an 

observation  vector,  {u}.  If  r and  u denote  prior  (initial)  estimates  of 

P P 

the  elements  of  these  vectors  based  upon  analytical  models,  then  the  "true” 
values  r and  u are  assumed  to  be  related  through  a sensitivity  matrix  by 
the  equation  (linearized  Taylor's  series  expansion) 

4-2 


It  is  important  to  note  at  this  point  that  the  sensitivity  matrix 
[T]  = [3u/3r]  and  the  estimator  matrix  [G]  are  updated  as  the  computational 
algorithm  iterates.  The  Kalman  Filter  equations  are  recursive  in  nature, 
and  eventually  the  estimator  matrix  [G]  will  reach  a steady-state  solution 
(cf.  Reference  10).  This  approach  (using  a linearized  Taylor's  series) 
is  conceptually  similar  to  applying  Newton's  method  for  solving  nonlinear 
simultaneous  equations. 

Furthermore,  as  indicated  previously,  the  method  of  satistical 
parameter  estimation  applies  when  there  are  "noisy"  data.  Experimental 
error  is  represented  by  { e } in  equation  (4-3)  and  is  accounted  for  by 
the  covariance  matrix  [ S _ c ] of  equation  (4-6)  . The  reader  interested 
in  further  details  of  Statistical  Parameter  Estimation  may  wish  to  consult 
References  7 and  8. 

Previous  structural  application  (Reference  9 ) of  statistical 
parameter  estimation  used  natural  frequencies  and  mode  shapes  (i.e.,  modal 
deflections  at  particular  locations)  as  elements  of  the  observation  vector. 
Reference  11  uses  frequency-response  data  and  assumes  steady-state  harmonic 
excitation.  The  estimation  procedure  can  be  applied  in  either  the  time- 
domain  (cf.  Reference  12)  or  the  frequency  domain.  The  present  study  uses 
the  time-domain  approach,  since  the  mathematical  model  is  non-linear. 

4.3  BAYESIAN  ESTIMATION  EQUATIONS 

The  J.H.  Wiggins  Company  has  made  use  of  "Statistical  Parameter 
Estimation"  previously  (e.g.,  References  9 and  11),  and  developments  in 
this  area  have  continued  apace.  Thus,  Isenberg  (Reference  13)  recently 
prepared  a systematic  account  of  statistical  parameter  estimation,  least- 
squares,  weighted  residual  methods,  etc.  From  this  systematic  study,  there 
resulted  a set  of  Bayesian  estimation  equations,  which  are  presented  in 
this  section  of  the  report.  (These  results  have  been  abstracted  from 
Reference  13).  Details  of  Isenberg's  results  are  included  herein  as 
Appendix  A) . 

In  Bayesian  estimation,  one  is  given  a prior  estimate  of  the  para- 
meters r , along  with  the  associated  covariance  matrix  S . One  then 
o ^ rr 

seeks  to  minimize  the  objective  function 
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where 


n 

£ 

n 

£ 

i=l 

j = i 

n 

n 

£ 

£ Vj 

II 

h-* 

j=i 

w = S 


-1 

ce 


(4-7) 


and 


w = S 


-1 

rr 


and  Srr  is  as  symmetric  matrix  as  is  S£e. 

The  variables  un  , u- , w. etc.,  are  defined  in  Appendix  A.  To  ori- 
ui  1 1 J 

ient  the  reader,  note  that  (u0^_ui)  represent  the  difference  between  experi- 
ment (u0£>serve<3)  and  the  theory  (u,.,  of  the  model).  The  terms  w^  repre- 
sents a "weighting"  matrix,  where  some  data  points  are  weighted  more  heavily 
than  others.  The  reason  for  this  "weighting"  is  that  all  the  data  points 
may  not  have  been  obtained  with  equal  certainty  (i.e.,  the  experimenter  may 
have  more  confidence  in  some  data  points  than  in  others) . For  example,  the 

experimenter  might  find  it  easier  to  measure  the  maximum  deflection  (6  , 

max 

say)  that  the  deflection  at  a (finite)  early  time  (e.g.,  at  t= . 1 millisec) . 
The  introduction  of  w^^  allows  the  experimenter  to  influence  the  model 
(which  is  being  estimated)  by  selectivity  weighting  his  data  points. 


Similarly,  equation  (4-7)  contains  the  parameters  (r  - r ■)  , where 

t h ^ i 

r_  is  the  initial  (a  priori)  estimate  of  the  i—  parameter.  The  weighting 
iirix,  • allows  the  structural  analyst  to  influence  the  parameters  esti- 
mated by  the  algorithm  in  a fashion  similar  to  the  experimentalist.  Thus, 
if  the  analyst  has  a great  deal  of  confidence  in  the  i—  parameter,  he  can 


o 

ma 


express  that  confidence  through  the  weighting  matrix,  w 
tice,  the  matrices  w and  w are  given  by 
.-1 


*-5 


w = S 


ee 


In  actual  prac- 


( 4 - 8a ) 


ana 


w = S 


-1 


rr 


( 4- 8b ) 
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If.  . , , 


The  user 


where  S and  S are  "covariance"  matrices  defined  in  Appendix  A. 
cz  rr 

of  PEBLS  controls  the  "route"  the  program  takes  (to  minimize  F in  equation 

4-7)  by  his  selection  of  S and  S_  . 

J zz  rr 

To  summarize  the  technique  of  Bayesian  estimation,  we  iteratively 
solve  the  sequence 


T = T <r  ) 

S = (w  + TTwT)-1  = S 

G = STTw  = G (re) 
r = rQ  + G [uQ-ue-T 


(re> 


(r  -r  ) ] 
o e 1 


( 4 -9a ) 
( 4 — 9b ) 

(4-9c) 
( 4 — 9 d ) 


and  repeat  the  next  sequence  with  r replacing  r . We  continue  to  iterate 
until  the  difference  between  r and  rg  becomes  sufficiently  small  and  the  se- 
quence has  therefore  converged.  The  final  value  of  S is  S*r,  the  revised 
parameter  covariance  matrix.  It  should  be  noted  that  when  the  response  is 
linear  with  respect  to  the  parameters,  only  one  iteration  is  required. 


As  discussed  in  Appendix  A,  the  matrix  T is  a "sensitivity"  matrix, 
i.e.,  the  change  in  response  due  to  a unit  change  in  a parameter.  (See 
equation  (42)  in  Appendix  A.)  The  matrix  S is  a covariance  matrix,  and  G 
is  the  "estimator"  matrix.  The  vector,  r,  is  the  set  of  parameters  (e.g., 
oc 0 , a^,  a2,  etc.)  of  the  mathematical  model.  (See  Section  3.0)  Equations 
(?-9)  are  used  within  PEBLS  to  obtain  the  revised  estimates  of  the  model 
parameters . 
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5.0 


RESULTS  FOR  THREE  (WIDELY  DIFFERENT)  TYPES 
OF  REINFORCED  CONCRETE  SLABS 


5.1  INTRODUCTION 

Static  and  dynamic  data  from  three  series  of  tests  were  used  in 
this  study.  The  results  of  the  modeling  are  contained  in  this  section  of 
the  report,  with  Section  5.2  devoted  to  Watt's  (deep)  slabs  (Reference  3), 
Section  5.3  to  Brown  and  Black's  conventional  slabs  (Reference  5),  and 
Section  5.4  on  Keenan's  laced  slabs  (Reference  6).  The  sensitivity  of 
the  results  is  given  in  Section  5.5,  and  additional  results  are  given  for  a 
three-parameter  model  (of  Watt’s  Slabs),  a four-parameter  model  (of  Brown 
and  Black’s  Slab),  and  a simple  model  of  Keenan's  Slabs. 

5.2  WATT'S  DEEP  SLABS  (Reference  3) 


Reference  3 reports  on  static  and  dynamic  tests  of  deep  R/C  slabs 
(span-to-thickness  ratio  of  4.12).  Two  types  of  slabs  were  tested  (namely 
mach-modeled  slabs  and  environmentally-modeled  slabs)  but  static  test  data 
were  available  only  for  the  latter.  The  slabs  were  square  in  plan,  but 
they  were  placed  over  a circular  cavity  and  thus  had  an  unsupported  diameter 
of  24  inches  in  the  center.  Dynamic  pressure  loads  were  applied  by  the  WES 
Small  Blast  Load  Generator,  which  involved  detonating  an  explosive  mixture 
of  gases  and  driving  a shock  wave  against  the  face  of  the  slab. 


The  dynamic  response  was  recorded  using  accelerometers  mounted  in 
the  center  of  the  slab  and  also  on  the  support  structure.  The  accelerometer 
traces  were  digitized  and  then  integrated  in  time  to  give  velocities  and 
displacements.  Pressure-time  histories  were  also  recorded.  For  further 
details,  see  Reference  3. 


The  computer  program,  PEBLS,  allows  velocities  to  be  part  of  the 
"observation  vector,"  {Au}.  Since  velocities  were  available  in  Watt's  test 
data,  they  were  used  (along  with  the  displacements)  to  form  the  observation 
vector.  Results  were  obtained  only  for  the  Environmentally  Designed  Slabs, 
test  specimens  EDI,  ED2  and  ED3.  Figure  5-1  shows  the  experimental  data 
which  was  used  as  input  to  the  program  and  Figure  5-2  gives  the  correspond- 
ing pressure-time  histories.  The  experimental  data  used  from  Figure  5-1  is 
tabulated  in  Table  I for  the  three  test  specimens,  EDI,  ED2  and  ED3 . 


In  addition  to  the  data  on  displacement  and  velocities,  note  that 
the  user  must  input  an  estimate  of  the  error  in  the  experimental  observa- 


tions. (The  error  in  observations  relates  to  the  matrix  [ S c c ] - see 


Equation  (4-8a)  discussed  previously).  Thus,  the  percent  errors  in  the 
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TABLE  I 


A.  EXPERIMENTAL  RESPONSE  FOR  THREE  SLABS 
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B.  CONFIDENCE  IN  THE  EXPERIMENTAL  OBSERVATIONS 
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velocities  and  displacements  are  also  given  in  Table  I (the  bottom  half) 
along  with  the  corresponding  variances.*  A small  percent  error  indicates 
good  confidence  in  the  data,  and  a large  error  indicates  poor  confidence. 
Thus  Table  I shows  that  we  have  more  confidence  in  the  displacements  at  t^ 
through  t7  (the  later  times)  than  we  do  in  those  near  the  initial  times 
(t^  through  t-j)  . 


The  results  of  the  parameter  estimation  calculations  are  shown  in 


m 


Table  II,  which  gives  the  ten  parameters  which  were  estimated  (P^^,  d^ 

Eplast'  etc • ) ' showing  their  initial  values,  their  values  estimated  for  each 

test  slab,  etc.  For  example,  the  initial  ("prior")  model  used  a yield  stress 
2 

of  525  LB/IN  (P, . ),  and  after  the  estimation  procedure  (for  slab  EDI)  the 
X in  2 

estimated  value  of  the  yield  stress  was  450  LB/IN  . Perhaps  the  more  impor- 
tant thing  to  note  in  Table  II,  however,  is  that  the  RMS  deviation  decreased 
from  .0428  initially  (Slab  EDI)  to  .0048  (when  the  revised  parameters  were 
used),  i.e.,  a factor  of  ten  improvement  in  the  RMS  value. 


The  deflection  vs.  time  is  shown  plotted  in  Figure  5-3  for  Slab  EDI, 
where  the  three  curves  represent 

(i)  the  experimental  data, 

(ii)  the  initial  model  response, 

and,  (iii)  the  model  response  using  the  "optimum"  parameters  just  esti- 
mated. 


TIME,  MILLI -SECONDS 

71-1  799 


Figure  5-3.  Deflection  - Time  History  of  Slab  ED  1, 

Showing  the  Initial  and  Final  Model  Results. 


♦The  variance  is  equal  to  the  square  of  the  standard  deviation.  See 
Reference  14. 
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It  is  also  noteworthy  that  the  static  load-deflection  curve 

(estimated  by  this  procedure)  is  a fairly  good  fit  to  the  static  test 

data.  (See  Figure  5-4).  These  results,  presented  in  Table  II  and  Figure 

5-3,  were  based  upon  a relatively  "good"  prior  model.  That  is,  given  the 

experimental  static  load-def lection  curve  of  Figure  5-4,  one  could  make  a 

pretty  fair  guess  about  the  values  of  the  parameters  P. . , d , E . , 

l in  nn  pisst 

etc . 

The  question  naturally  arises,  "What  if  one  didn't  have  the  static 
load-deflection  curve?"  In  this  case,  he  would  rely  on  a combination  of 
experience  and  analysis  to  come  up  with  a prior  model  for  the  non-linear 
load-deflection  characteristic.  But  then  the  question  is,  "How  well  would 
the  estimation  procedure  perform,"  or  more  specifically,  "Will  the  estima- 
tion procedure  still  arrive  at  a good  bi-linear  approximation  to  the  exper- 
imental load-def lection  curve?" 

To  test  this  problem  out,  a poor  initial  model  (with  a yield  stress 
2 

of  only  200  LB/IN  ) was  input  to  the  program  along  with  the  experimental 
response  data  from  Table  I.  The  results  are  shown  in  Table  III  and  Figures 
5-5  and  5-6.  Referring  to  Figure  5-5,  note  that  the  response  curve  for  the 
revised  model  matches  the  experimental  data  almost  exactly,  although  the 
initial  (prior)  model  gives  rather  poor  results. 

Figure  5-6  shows  similar  results,  where  the  load-deflection  curve 
fox  the  final  model  is  much  closer  to  the  experimental  data  than  is  the 
initial  model.  These  preliminary  results  were  encouraging,  and  it  appeared 
that  the  1 d-o-f  model  (as  used  herein)  was  able  to  match  the  experimental 
response  for  deep  slabs. 

Refer  again  for  a moment  to  Table  II.  Note  that  Table  II  contains 
both  the  initial  estimates  of  error  in  the  parameters  as  well  as  revised 
error  estimates.  For  some  parameters  (e.g.,  EMUD) , the  initial  error 
estimate  was  50%,  and  after  data  for  three  slabs  had  been  supplied  to  the 
computer  program,  the  error  was  barely  reduced  (i.e.,  the  revised  error  was 
49%).  However,  note  that  other  parameters  (like  PLIN,  for  example)  began 
with  an  estimated  error  of  50%,  which  was  subsequently  reduced  to  17%, 
i.e.,  a significant  improvement. 

Thus,  it  seemed  that  everything  was  working  exactly  as  had  been 
intended,  and  the  parameters  (see  Table  II)  were  coming  out  with  "reasonable 
values.  At  the  suggestion  of  the  Technical  Monitor,  we  simplified  the 
model  from  ten  parameters  (Table  II)  to  four  parameters,  namely: 
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Figure  5-4.  Pressure  Versus  Midspan  Deflection:  Test  Data  and  Revised  Model 


TABLE  III 


A.  PARAMETER  VALUES,  BEGINNING  WITH  A POOR  MODEL 


ECTION,  INCHES 


Initial  Guess,  and  Final  Estimate 


• A single  mass  parameter,  u 

avg 

• Two  spring-force  parameters  (P, . and  d ) 

lin  lin 

• /.nd  one  force  parameter,  pscale  (the  a value  was  set 
equal  to  unity) . 

The  results  of  this  "four-parameter  model"  are  shown  in  Table  IV, 
and  Figures  5-7  through  5-9,  respectively,  for  slabs  EDI,  ED2,  and  ED3. 

Note  that  by  adjusting  these  four  parameters  (p,  P-^n>  d}.in  and  ^scale^' 

The  parameter  estimation  program  was  still  able  to  closely  fit  the  dis- 
placement time-histories.  The  parameters  still  seemed  to  have  reasonable 
values,  but  it  was  somewhat  disturbing  that  the  force-deflection  para- 
meters (p, . , d.,  . ) did  not  match  the  static  load-deflection  data  as  well 
as  had  been  expected.  These  discrepancies  were  thought  to  be  due  to  slight 
variations  in  the  individual  slabs,  however. 

At  this  point  in  the  study,  everything  looked  rosy,  and  it  was 
(somewhat  naively)  assumed  that  similar  results  would  be  obtained  for  the 
other  types  of  slabs  as  well.  However,  that  was  not  the  case,  as  is  dis- 
cussed in  the  sections  which  follow. 

5.3  RESULTS  FOR  BROWN  AND  BLACK'S  CONVENTIONAL  R/C  SLABS  (Reference  5) 

Reference  5 provides  both  static  and  dynamic  test  data  on  conven- 
tional R/C  slabs,  where  the  slabs  were  square  in  plan  and  had  a span  of  29 
inches  and  a thickness  of  .89  inches.  The  static  test  data  are  shown  in 
Figure  5-10.  Note  that  these  static  tests  were  continued  to  destruction  of 
the  slab,  i.e.,  until  the  slab  could  carry  no  more  load,  at  a deflection  of 
approximately  4 inches.  At  this  point  in  the  test,  the  concrete  was  highly 
fractured,  and  the  reinforcing  steel  was  stretched  to  its  ultimate  capac- 
ity. (For  additional  details  and  photograph  of  the  damaged  slabs,  see  Re- 
ference 5)  . 

Brown  and  Black's  tests  are  noteworthy  for  the  fact  that  they 
determined  the  ultimate  pressure  capacity  of  their  slabs  (From  20  to  25 
psi,  see  Figure  5-10)  and  the  corresponding  maximum  deflection  (i.e., about 
4 inches).  This  information  was  not  obtained  for  Watt ' s deep  slabls,  dis- 
cussed previously  (cf.  Figure  5-6).  Thus,  from  their  static  tests,  Brown 
and  Black  knew  that  a deflection  of  about  4 inches  corresponded  to  the 
ultimate  capacity  of  their  slabs.  Yet,  in  the  dynamic  tests  (see  Figure 
5-11),  the  slabs  were  all  driven  to  amplitudes  well  in  excess  of  5 inches. 

Not  surprisingly,  each  of  the  dynamically  tested  slabs  was  severely 
damaged.  (See  Reference  5 for  photos  of  the  tested  slabs).  These  tests 

5-11 


Table  IV 

Estimated  Parameters  for  Watt's  Slabs 


PARAMETER 


MASS  COEFFICIENT  y 

.93 

1.00 

.95 

YIELD  FORCE  P1jf} 

400 

360 

443 

ELASTIC  DISPLACEMENT  dlin 

.137 

.173 

.168 

PRESSURE  SCALE  FACTOR  Pscale 

1 

.96 

1 

.82 

.90 

RMS  DEVIATION  ARMS 


005 


.015 


.003 


Deflection  Time-History  Fitted  by  Four-Parameter  Model. 
(Watt's  Slab  ED-1,  Reference  3) 


INITIAL  ESTIMATE 


Deflection  Time-History  Fitted  by  Four-Parameter  Model. 
(Watt's  Slab  ED-2,  Reference  3) 


(by  Brown  and  Black)  were  therefore  fundamentally  different  than  Watt's 
tests,  where  the  ultimate  capacity  of  the  slabs  was  not  measured,  and  the 
dynamic  tests  were  run  at  less  than  the  ultimate  strength.  On  the  other  hand, 
the  authors  of  Reference  5 did  measure  the  ultimate  strength  of  their  slabs, 
and  then  they  tested  (dynamically)  to  well-beyond  the  static  capability. 

These  basic  differences  in  the  test  philosophy  and  test  results  were  not 
immediately  recognized  by  the  writer.  It  was  only  after  difficulties  were 
experienced  with  the  parameter  estimation  procedure  that  a closer  look  was 
made  into  the  test  details. 


Initially,  a 10-parameter  model  was  used  to  fit  the  test  data  for 

slab  ID-2.  (The  pressure-time  history  for  this  test  is  given  in  Figure 

5-12,  and  the  displacement  time-history  is  shown  in  Figure  5-13.)  This 

multi-parameter  model  resulted  in  a good  fit  to  the  time-history  data  (see 

Figure  5-14)  and  the  force-deflection  results  also  looked  promising  (cf. 

Figure  5-15).  However,  in  contrast  with  the  earlier  results  (Section  5.2) 

the  parameters  for  a,  p , etc,  did  not  appear  to  be  physically  realistic. 

To  get  a better  understanding  of  what  was  happening,  the  10-parameter  model 

was  simplified  (e.g.,  p was  held  constant,  a was  held  constant,  etc.)  but 

E . (the  post-yield  slope)  was  still  allowed  to  vary  (since  it  was 
piast 

expected  that  a bi-linear  characteristic  would  be  needed  to  approximate  the 
force-deflection  curve) . Several  parameter  estimation  runs  were  then  made 
(they  are  tabulated  in  Table  V)  each  with  the  same  result:  the  program 
arrived  at  a "nearly-linear"  force-deflection  curve  when  fitting  the  dynamic 
test  data.  The  fit  to  the  dynamic  data  was  fairly  good  in  each  case,  as 
denoted  by 


A 

rms 


N 

I X 

1=1 


1/2 


which  is  a measure  of  the  deviation  between  the  model  (theory)  and  the  test 
(experiment) . 

At  this  point,  it  was  decided  to  try  to  use  a linear  force-deflection 
curve  and  see  what  the  parameter  estimation  program  would  do.  Again,  it 
achieved  a fairly  good  fit  to  the  dynamic  test  data,  with  Arms  less  than  .01. 
(See  the  last  line  of  Table  V.)  These  results  were  disturbing,  at  least 
initially,  since  they  were  unexpected.  What  had  been  expected  was  that  some 
form  of  strain-hardening  (e.g.,  like  Figure  5-15)  would  be  required  to 
satisfactorily  fit  the  dynamic  data.  It  was  thought  that  the  strain- 
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Figure  5-15.  Multi-Parameter  Model,  Showing  Fit  to  the  Static  Force 
Deflection  Data 


Stic  slope 


urdening  would  be  required  because  the  static  load-deflection  data  was 
strongly  non-linear  (c.f.  Figure  5-10).  However,  the  parameter  estimation 
equations  kept  fitting  the  dynamic  data  and  using  a nearly-linear  model  to 
do  so . 


Final ly,  after  some  hand  calculations  using  the  method  of  equivalent 
linearization  (Reference  15)  it  was  recognized  that  the  dynamic  data  (or  at 
least  the  maximum  deflection,  6max)  can  always  be  fitted  by  an  "equivalent 
linear  model".  For  any  loading  function,  P(t),  the  response  of  the  1 d-o-f 
system  is  well-approximated  by 


max 


is  the  frequency  of  the  equivalent  linear  oscillator 


and  k is  the 
eq 


If  one  begins  with  the  non-linear  differential  equation 


and  the  approximate  sdution 


max 


(where  f (y) is  any  non-linear  force-deflection  curve) then  multiplying  both 
nl  — 

sides  of  ehe  equation  (5-3)  by 


y = w 6 sin  w t 
1 eq  max  eq 

and  integrating  the  result  gives 


eq  max 


The  limits  on  the  integration  are  from  y=0  to  y=6  , or  in  terms 

r max 

of  the  time,  t,  from  t=0  to  t = . At  both  these  times,  the  velocity,  y 

is  zero,  and  the  kinetic  energy  term  vanishes.  Equation  (5-5)  then 


eo  max 


which  says  that  the  strain  energy  (stored  in  the  nonlinear  spring)  is  equal 
to  the  "work  done"  by  the  applied  force  p(t). 

As  a first  approximation,  one  might  try  p(t)  = p^  = a constant,  if 
the  loading  function  does  not  change  much  during  the  initial  response.  Then, 
the  right-hand-side  of  equation  (5-6)  gives 

/P  w 6__  sin  w__t  dt  = 2p^w  6 (5-7) 

*0  eq  max  eq  ro  eq  max 

and  one  finally  has 
f max 

/ f„i(y)dy  = 2p  w 6 (5-8) 

J nlw'  1 ro  eq  max 

0 

which  can  be  solved  for  the  equivalent  linear  frequency  (w  ) and  thence 

the  equivalent  linear  stiffness  (''eg)  • Note  that  ueg  is  given  in  terms  of 

the  area  (strain  energy)  under  the  force-deflection  curve,  and  the  maximum 

displacement,  £max-  Also  note  that  no  restriction  has  been  p’aced  on  the 

force-deflection  characteristic,  f (y)  . 

nl 

Based  upon  the  foregoing  discussion,  it  is  the  author's  opinion  that 
the  parameter  estimation  algorithm  performed  an  "equivalent  linearization" 

(or  something  similar)  to  achieve  a good  fit  to  the  dynamic  test  data  of 
Reference  5 with  a nearly-linear  force-deflection  curve.  If  this  hypothesis 
is  correct,  it  is  only  natural  to  inquire  as  to  why  something  similar  did 
not  occur  previously,  e.g.,  in  the  case  of  Watt's  slabs  (Section  5.2).  A 
possible  explanation  for  this  difference  in  behavior  is  that  Watt's  slabs 
were  tested  in  a "moderately  nonlinear"  range  (with  the  maximum  deflection 
ax  equal  to  just  2 or  3 elastic  deflections,  d^n)  whereas  Brown  and 
Black's  slabs  were  loaded  to  their  ultimate  capacity.  For  the  latter,  the 
maximum  deflections  were  on  the  order  of  15  to  20  times  the  elastic  deflec- 
tions. As  stated  previously,  there  was  a fundamental  difference  in  the 
testing  philosophies  used  in  References  3 and  5.  This  point  will  be  dif- 
cussed  in  more  detail  in  Section  5.9. 

5.4  RESULTS  FOR  KEENAN'S  "LACED"  R/C  SLABS  (REFERENCE  6) 

Reference  6 presents  results  for  static  and  dynamic  tests  of  uni- 
formly-loaded one-way  reinforced,  "laced"  slabs.  The  slabs  had  a span  of 
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ypical  Longitudinal  and  Transverse  Section  of  Laced 
lab  {Keenan,  Reference  6) 
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Figure  5-19.  Variation  of  Deflection  Across  the  Semi-Span 
(Keenan's  Slab  D3-1) 
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6 feet  (72  inches)  and  a thickness  of  6 inches,  giving  a span- to-aepth 
ratio  of  12.  These  slabs  are  unusual  because  of  their  "lacing"  reinforce- 
ment, which  is  shown  in  Figure  5-16.  Keenan’s  tests  are  also  unusual  by 
virtue  of  the  fact  that  the  tests  were  conducted  at  several  different 
pressure  levels.  Recall  that  Watt's  tests  (Reference  3)  and  Brown  and 
Black's  (Reference  5)  were  each  conducted  at  basically  a single  (dynamic) 
pressure  level.  Keenan  recognized  that  a nonlinear  system  (such  as  an  R/C 
slab)  should  be  tested  throughout  its  nonlinear  range.  Thus,  his  tests 
began  at  a low  level  of  (dynamic)  pressure,  p(t),  and  became  progressively 
more  severe.  However,  a testing  program  of  this  type  requires  either  several 
identical  specimens  (which  is  expensive)  or  several  tests  on  the  same 
specimen  (which  results  in  cyclic  loading  and  progressive  failure) . Both 
alternatives  are  undesirable,  from  the  point  of  view  of  the  test  engineer. 
Keenan  chose  the  latter  course,  and  he  successively  tested  the  slabs  to 
higher  and  higher  pressures.  The  most  significant  data,  from  the  point  of 
view  of  the  present  study,  occurred  on  the  "first  shot."  Thus,  the  atten- 
tion is  focused  herein  on  Keenan's  test  specimen  D3-1.  The  pressure  vs. 
time  history  and  displacement  time  history  for  this  slab  are  shown  in 
Figure  5-17.  A typical  static  pressure  vs.  deflection  curve  for  such  slabs 
is  shown  in  Figure  5-18. 

Referring  to  Figure  5-17,  note  that  there  are  three  displacement 
traces,  namely  Dl,  D2,  and  DC.  These  traces  refer  to  three  separate  dis- 
placement-measuring gages,  mounted  at  points  along  the  span  of  the  slab, 
with  gage,  DC  in  the  center.  By  using  three  displacement  gages  (instead  of 
just  one,  in  the  center)  Keenan  observed  that  the  displacements  varied  in 
space , as  well  as  with  time.  For  example,  Figure  5-19  shows  the  readings  of 
the  three  displacement  gages  (connected  by  straight  lines)  at  identical 

It  is  clear  from  Figure  5-19  that  more  than  one  mode  is  participa- 
ting in  the  response,  and  in  particular,  it  appear  that  the  first  mode 
(like  sin  — ) is  interacting  with  the  third  mode  (like  sin  — = — ) . 

Lt 

Several  attempts  were  made  to  "fit"  the  central  response  (gage  DC) 
of  Keenan's  slab  D3-1  to  the  one  degree-of-freedom  model.  The  results 
which  were  successful  are  given  in  Table  VI  in  terms  of  the  model  of 
parameters.  Of  particular  note  was  the  fact  that  the  "force  coefficient", 

. i , was  always  estimated  as  befinning  at  a very  low  value  (e.g.,  .02,  etc.) 

initially.  Attempts  to  fit  the  response  curve  with  a constant  a (like  an 
aavg  suc^  as  was  successful  with  Watt's  slabs)  were  totally  unsuccessful 
in  this  case. 
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By  beginning  with  a - .C2  and  then  going  up  to  = 1,  the  force 
coefficient  aeemphasized  the  initial  pressure  peak,  since  the  force  applied 
to  the  1 d-o-f  model  is  given  by  the  product  of  a and  p(t) . Without  this 
significant  "re-shaping"  of  the  pressure  pulse,  it  was  not  possible  to 
achieve  a very  good  fit  to  the  dynamic  response  data  for  Keenan’s  laced 
slabs . 

Keenan's  tests  used  prima-cord  explosives  placed  in  vented  firing 
tubes  to  produce  the  dynamic  over-pressure  on  the  slabs.  Possibly,  because 
of  non-uniformity  of  pressure  loading  or  perhaps  because  of  gas  leakage 
(around  the  edges  of  the  slab)  Keenan's  initial  pressure  reading  may  have 
been  too  high.  Or,  another  possibility  is  that  the  interaction  of  the 
fundamental  mode  with  the  higher  modes  occurred  in  such  a fashion  that 
the  "effective  force"  on  the  fundamental  mode  was  initially  reduced. 

In  that  regard,  it  is  worth  noting  that  tha  slab  failure  occurred 
(both  statically  and  dynamically)  with  a hinging  action  at  the  supported 
ends  and  in  the  center  (see  Figure  5-20).  Yet  from  the  linear  theory  of 
vibrations,  one  expects  the  slab  to  respond  initially  (i.e.,  at  small 
amplitudes)  in  terms  of  the  vibration  modes  of  a clamped-clamped  beam. 

Thus,  the  slab  response  must  make  a (fairly  significant)  transition  from 
several  vibration  modes  at  small  amplitudes  to  a hinged  mode  at  large 
amplitudes.  This  modal  transition  might  be  the  "cause"  of  the  reduced  force 
coefficient  (a)  wiiich  the  estimation  algorithm  seems  to  require  in  fitting 
the  response  data  for  Keenan's  slab  D3-1. 
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SENSITIVITY  OF  THE  RESULTS  TO  MODEL  PARAMETERS 


In  the  parameter  estimation  algorithm,  the  "sensitivity  matrix", 

[T] , is  calculated  and  used  in  arriving  at  the  revised  parameters.  The 
individual  elements  of  [T]  constitute  the  "change  in  response  per  unit 
change  in  a given  parameter,"  and  these  sensitivities  are  tabulated  in 
the  computer  printout.  To  give  the  reader  a feeling  for  which  parameters 
have  the  most  significance,  Tables  VII  through  IX  have  been  prepared . These 
tables  give  the  sensitivities  at  three  locations  in  the  dynamic  response 
time-history,  namely 

(i)  the  sensitivities  during  the  initial  response,  when  the 
deflection  is  near  zero, 

(ii)  the  sensitivities  at  the  mid-response,  when  the  velocity  is 
a maximum,  and 

I (iii)  the  sensitivities  at  the  maximum  response. 

Admittedly,  Tables  VII  through  IX  are  somewhat  misleading,  in  that 

the  absolute  sensitivities (which  are  used  and  output  by  the  computer  program) 

-2 

have  been  tabulated.  Thus,  Table  VII  gives  Ay. /Ay  as  -2.4x10  , which  means 

“2 

the  displacement  y^  changes  1 2.4x10  j for  a unit  change  in  the  mass  coef- 
ficient y.  Note,  however,  that  y itself  may  have  been  only  1x10  initially 
; so  Ay^/Ay  could  mean  a relative  change  in  y^  of  240%.  It  would  be  desirable 

(if  possible)  to  present  the  relative  sensitivities,  like  (Ay^/y^)  (divided 
by)  (Ay/y) . To  do  so  at  this  point  in  time  would  require  many  hours  of 
searching  through  old  computer  runs  and  was  not  thought  to  be  worthwhile. 
Subsequent  changes  to  PEBLS  (to  print  out  relative  as  well  as  absolute 
sensitivities)  are  definitely  indicated  at  a future  date. 

This  is  not  to  imply  that  Tables  VII  through  IX  are  not  useful; 
note,  for  example  in  Table  IX  that  the  sensitivity  Ay^/Aao  was  75.5  at  the 
maximum  response  point.  Clearly  y^  is  very  sensitive  to  aQ,  as  indicated 
by  the  high  value,  namely  75.5. 

5.6  THREE-PARAMETER  MODEL  RESULTS  FOR  WATT'S  SLABS 

I As  stated  in  Section  2.4,  it  eventually  became  evident  that  a three- 

paiameter  model  (involving  X,  the  load-mass  factor,  P. . and  d.. , ) was 

jl  in  i in 

, the  simplest  elasto-plastic  model  one  can  consider  (cf.  Biggs,  Reference  2). 

Accordingly,  an  attempt  was  made  to  "fit"  Watt's  dynamic  test  data  using 
this  three-parameter  model.  The  results  are  summarized  in  Table  X,  which 
shows  that  a pretty  fair  fit  (e.g.,  an  RMS  error  of  .01)  can  be  obtained 
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Table  VII 

Response  to  Various  Model  Parameters  (Watt's  Slab  ED-1)* 


Table  VIII 


Typical  Sensitivities  of  the  Response  to  Various  Model  Parameters 
(Brown  and  Black's  Slab  ID-2)* 


INITIAL  RESPONSE 

(NEAR-ZERO 

DISPLACEMENT) 

MID-RESPONSE 
(NEAR  MAXIMUM 
VELOCITY) 

MAXIMUM  RESPONSE 
(AT  MAXIMUM 
DISPLACEMENT) 

SENSITIVITY  TO 

FORCE  PARAMETERS: 

,16 

6.5 

3.8 

Ay-j/AQj 

-1.5  x 10"4 

-.43 

-.41 

SENSITIVITY  TO 

MASS  PARAMETERS: 

AVApo 

1 x 10'3 

-1.0 

-.53 

Ay  -j/Ayd 

.16 

.34 

.19 

SENSITIVITY  TO 

BI -LI NEAR  SPRING 
PARAMETERS: 

AV4Plin 

“ ^ 
-2x10 

-.13 

-.30 

AVAdlin 

l.Ox  10-3 

.53 

1.08 

^i/AEpiast 

0 

-.06 

-1,05 

SENSITIVITY  TO 

PRESSURE-SCALING 

PARAMETER: 

AVAPscale 

.22 

2.5 

5.04 

r *See  also  Table  VI,  row  two. 

D * 
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Table  IX 


Typical  Sensitivities  of  the  Response  to  Various 
Model  Parameters  (Keenan's  Slab  D3-1)* 


INITIAL  RESPONSE 

(NEAR-ZERO 

DISPLACEMENT) 

MID-RESPONSE 
(NEAR  MAXIMUM 
VELOCITY) 

MAXIMUM  RESPONSE 
(AT  MAXIMUM 
DISPLACEMENT) 

SENSITIVITY  TO 

FORCE  PARAMETERS : 

AVAao 

.49 

27.3 

75.5 

Ay./Aad 

-.01  v 

-5.3 

-16.0 

SENSITIVITY  TO 

MASS  PARAMETERS: 

AVAuo 

-.013 

-1.46 

-4.36 

AVAud 

3.9  x 10"6 

.018 

.062 

SENSITIVITY  TO 
BI-LINEAR  SPRING 
PARAMETERS: 

Ayi/APlin 

-8.4 x 10"6 

-7.5  x 10”  3 

-4.5x  10"2 

Wi^lln 

4.3x  10"3 

2.67 

7.82 

4V4Eplast 

0 

-2.2  x 10-4 

-1.5  x 10"2 

SENSITIVITY  TO 

PRESSURE-SCALING 

PARAMETER 

«V4Pscale 

5.6  x IQ'3 

.95 

3.74 

*See  also  Table  V,  row  six. 
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with  the  simple  model.  Note,  also,  the  "yield  force"  was  typically 

estimated  at  about  450  Lb/in2,  in  pretty  fair  agreement  with  the  static 
test  data,  (cf.  Figure  5-4  and  5-21).  As  shown  in  Figure  5-21,  the  estimated 
elasto-plastic  load-deflection  curves  are  quite  close  to  the  results  from 
static  tests.  Thus,  one  can  get  a fairly  good  estimate  of  the  analytical 
load-deflection  curve  from  static  test  data. 

Note  that  Table  X also  gives  estimated  load-mass  factor  of  close  to 
unity.  A load-mass  factor  of  exactly  unity  is  predicted  if  the  slab  "punches 
through"  with  no  spanwise  bending.  Watt's  tests  (Reference  3)  showed  that 
many  of  his  slabs  did  "punch  through" , which  accounts  for  the  load-mass 
factors  estimated  and  shown  in  Table  X.  It  is  felt  that  the  parameters  given 
in  Table  X are  "physically  realistic"  for  Watt's  slabs,  and  this  confirms 
the  earlier  result  that  simple  models  will  work  pretty  well  for  Watt's  data. 


Table  X.  Three-Parameter  Model  Estimates 
for  Watt's  Slabs 


PARAMETER 

SLAB  EDI 

SLAB  ED2 

SLAB  ED3 

LOAD-MASS  FACTOR, 

00 

1.19 

1,05 

YIELD  FORCE, 

422 

439 

480 

ELASTIC  DISPLACEMENT, 

d.  . 

1 in 

CO 

.137 

.135 

GOODNESS  OF  TIT, 

i.e.,  Arms 

.006 

.016 

005 
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Pressure  Versus  Midspan  Deflection  Showing  Tes 
Range  of  Final  Estimates  Found  Using  PEBLS 


! r 
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FOUR-PARAMETER  MODEL  FOR  BROWN  AND  BLACK’S  SLAB  ID-2 


In  keeping  with  the  theme  of  using  the  "simplest  possible  model", 
a four-parameter  model  was  attempted  with  Brown  and  Black's  test  data.  The 


four  parameters  were  A (the  load-mass  factor)  P. . , d. . and  E . , , the 

lin  lin  plast 


latter  being  added  primarily  because  of  the  stronger  non-linear  force- 
deflection  results  (cf.  Figure  5-10). 


Table  XI. 

Results  of 

Estimation 

Procedure 

for  Brown  ai 

COMPUTED  AND 

OBSERVED  DISPLACEMENT  RESPONSE 

POINT 

TEST  TIME 

U-ANALTT1C 

U-OBSERVE  D 

I*U-UTEST 

1 

0 . 2500  00*0  2 

0.94 J41E-01 

0. BOOOOE-Ol 

0.14341 E-0 1 

2 

0.50000E-02 

0.54926 

0.50000 

0.49256E-01 

3 

0.75000E-02 

1.5285 

1 .2000 

0.32851 

4 

0, 10000E-01 

2.9738 

2.0000 

0.97379 

5 

0.12500E-01 

4,6882 

3 • 00O0 

1.6882 

0 

0. 15000E-01 

6.4588 

4.2000 

2.2588 

7 

0.17500E-01 

8.0667 

5.4000 

2.o667 

8 

0.20000E-01 

9.2944 

6. 1000 

3.1944 

9 

0.22500E-01 

9.9598 

6.8000 

3.1598 

10 

0. J0000E-01 

9.6021 

7.0400 

2.5621 

OBJECTIVE  DUE  TO  OBSERVATIONS 
OBJECTIVE  DUE  TO  PARAMETERS 


179.769 

0.000000 


RMS  ERROR  IN  OBSERVATIONS 
TOTAL  OBJECTIVE  EUNCTION 


2.07064 
1 7*. 7e9 


COMPUTED  AND  OBSERVED  DISPLACEMENT  RESPONSE 


>0INT 

TEST  TIME 

ii-analitic 

U-0BSERVED 

Y«U-UTEST 

1 

0.25000E-02 

0.80B62E-01 

0. 80000 E 

-01 

0.86156t-03 

2 

0.50000E-02 

0.47053 

0.50000 

-0. 29471E-01 

3 

0.75000E-02 

1.1535 

1.2000 

-0.46527E-01 

4 

0. 10000 K “01 

2.OB01 

2.0000 

0.801 2UE-01 

5 

0. 12500E-01 

3.1507 

3.0000 

0. 15069 

6 

0. 15000E-01 

4.2646 

4.2000 

0.64625E-01 

7 

0. 17500F-01 

5.3191 

5.4000 

-0.80906E-01 

e 

0. 20000 E “01 

6 .2040 

b.  1000 

0.10399 

9 

0.22500E-01 

6.8204 

6.0000 

0 • 20403E-01 

10 

0.30000E-01 

6.8452 

7.0400 

-0.19475 

OBJECTIVE.  DUE  TO  OBSERVATIONS  s 

1 .06444 

RNS 

error  in  observations 

OBJECTIVE  DUE  TO  PARAMETERS 


1 .53308 


TOTAL  OBJECTIVE  EUNCTION 


0.959385E-01 

2.59753 


Table  XTI.  Parameters  Estimated  for  Brown  and  Black's  Slab 


PARAMETER 

confidence  LEVELS 

INITIAL 

LATEST 

PERCENT 

PARAMETER  ESTIMATE 

ESTIMATE 

VARIANCE 

CON* IDENCE 

PLIN 

14.000 

19.266 

2.8877 

12.138 

DLIN 

0.25000 

0.33194 

0.1 1518E-01 

42.928 

EPLS 

1.4000 

0.37343 

0.33278 

41.205 

EMU1 

0.70000 

0.82844 

0.36831E-01 

27.416 
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The  results  from  the  computer  run  are  shown  in  Table  XI  which 
gives  the  initial  response  (top  half)  and  the  final  response  (bottom  half 
of  Table  XI) . Note  that  the  RMS  error  in  observation  went  from  a value  of 

2.07  (initially)  to  a value  of  only  .0959  (final  estimate). 

The  final  estimated  parameter  values  are  given  in  Table  XII, 
which  shows  P^n<  ^iin'  Eplast  anc^  ^ (labeled  EMUl)  each  having  a 
significant  improvement  in  confidence  level  (e.g.,  from  75%  to  12%,  etc.) . 
Note  also  that  the  values  estimated  appear  to  be  "physically  reasonable"  and 
acceptable  from  an  intuitive  standpoint.  (For  example,  see  the  text  by  Biggs, 
Reference  2,  for  typical  load-mass  factors,  and  Figure  5-22  for  the  static 
load-deflection  curve).  Biggs  test,  p.  214,  gives  theoretical  load-mass 
factors  in  the  range 

.51  < X < .67 

for  square  two-way  slabs  (a/b  = 1).  Table  XII  gives  a load-mass  factor  of  .82, 
which  is  in  fair  agreement  with  the  theory.  Note  that  Biggs’  values  are 
calculated  assuming  "yield-line  theory"  for  the  failure  of  the  slab. 

Figure  5-22  shows  the  estimated  bi-linear  pressure  vs.  deflection 
fair  amount  of  scatter  is  indicated  by  the  test  data,  suggesting  that  the 
slabs  were  of  non-uniform  quality.  Again  one  sees  that  simplified  theory 
and  the  static  load-deflection  data  provide  fairly  good  estimates  of  the 
final  model  parameters. 

5.8  FOUR  PARAMETER  MODEL  FOR  KEENAN'S  SLAB  D3-1 

Keenan's  data  (on  his  Slab  D3-1)  proved  to  be  fairly  difficult  to 
"fit"  with  the  simple  four-parameter  model.  Several  computer  runs  were  made 
using  the  data  shown  in  Figure  5-17).  Discussed  previously  in  Section  5.4). 
Accepting  the  data  as  accurate,  the  computer  program  PEBLS  originally  esti- 
mated a load-mass  factor  of  about  three  (3)  for  Slab  D3-1.  This  result  was 
unexpected  (and  we  felt  it  was  incorrect)  so  a closer  look  was  taken  at  the 
data  (Figure  5-17). 

Referring  to  Figure  5-17,  note  that  tne  displacement  trace  DC  shows 
a deflection  of  nearly-zero  out  to  a time  of  .002  sec  (2  milli-seconds) . 

Note  also  that  the  pressure,  p ( t ) is  the  highest  during  this  early  time 
period.  Thus  the  possibility  occurred  that  perhaps  Keenan's  displacement 
trace  was  in  error,  since  when  the  pressure  is  greatest,  (and  the  accel- 
eration is  a maximum)  one  would  expect  that  displacement  would  occur  (and 
not  be  near-zero) . 
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Table  XIII 


Results  of  Estimation  Procedure  for  Keenan's  Slab  D3-1 


COMPUTED  AND 

OBSERVED  DISPLACEMENT  RESPONSE 

POINT 

TEST  TIME 

U-ANALIT1C 

U-OBSERVEO 

KaU-UTEST 

1 

0.20000E-02 

0.32253 

0.20000 

0.12253 

2 

0. 30000F-02 

0.60658 

0.35000 

0.25658 

3 

0.40000E-02 

0.89463 

0.53000 

0.36463 

4 

0.50000E-02 

1.1630 

0.75500 

0.40801 

5 

0.60000E-02 

1.3919 

0.96000 

0.43189 

6 

0.70000E-02 

1.5662 

1 .1600 

0.40623 

7 

0.80000E-02 

1.6745 

1.2900 

0.38446 

8 

0.10600E-01 

1.5771 

1.4600 

0.11710 

OBJECTIVE  DUE  TO  OBSERVATIONS  a 83.9786 
OBJECTIVE  DUE  TO  PARAMETERS  a 0.000000 


RMS  ERROR  IN  OBSERVATIONS  a 0.334168 
TOTAL  OBJECTIVE  FUNCTION  a 83.9786 


COMPUTED  AND  OBSERVED  DISPLACEMENT  RESPONSE 


POINT 

TEST  TIME 

U-ANALYTIC 

U-OBSERVED 

YsU-UTEST 

1 

0.20000E-02 

0.19966 

0.20000 

-0.34455E-03 

2 

0. 30000E-02 

0.37947 

0.35000 

0.29471E-01 

3 

0.40000E-02 

0.57767 

0.53000 

0.47666E-01 

4 

0.50000E-02 

0.77838 

0.75500 

0.2337BE-01 

5 

0.60000E-02 

0.96830 

0.96000 

0.83005E-02 

6 

0.70000E-02 

1 .1374 

1.1600 

-0.22642E-01 

7 

0.80000E-02 

1.2778 

1.2900 

-0.12229E-01 

8 

0.10600E-01 

1.4548 

1.4600 

-0.52067E-02 

objective  due  to  observations  a o. 459336 
OBJECTIVE  DUE  TO  PARAMETERS  a 1.98882 


RMS  ERROR  IN  OBSERVATIONS  a 0.236729E-01 
TOTAL  OBJECTIVE  FUNCTION  a 2.44815 


Table  XIV  . 


Parameters  Estimated  for 


Keenan's  Slab  D3-1 


parameter 

CONFIDENCE  levels 

INITIAL 

latest 

PERCENT 

parameter  estimate 

ESTIMATE 

variance 

CONFIDENCE 

PLIN 

110.00 

86.471 

112.11 

9.6256 

DLIN 

0.34000 

0.88527E-01 

0.23315E-01 

44.909 

EPLS 

1 .0000 

0.75941 

24.797 

497.97 

EMU! 

0,80000 

I.l6b3 

0.98243E-01 

39.180 

For  this  reason,  the  displacement  data  was  "shifted"  to  the  left 


in  time  a distance  of  .002  sec.  Thus,  the  time  t = 0 was  taken  as  the  time 
that  the  pressure-trace  began,  and  the  time  on  the  displacement  trace  was 
shifted  to  a new  time,  (t'  = t - .002).  When  this  adjustment  to  Keenan's 
data  was  made,  a fairly  satisfactory  "fit"  was  obtained  using  the  four- 
parameter  model.  * 

1 

Using  the  static  load-deflection  curve  (Figure  5-18) , one  can 

estimate  P. . , d, . , and  E . . fairly  well  for  Slab  D3-1.  Then,  using 

an  lin  plast  ^ 

these  values  with  a representative  load-mass  factor,  the  results  shown  in 

Table  XIII  are  obtained  (top  half).  Note  that  the  RMS  error  value  is  .33, 

indicating  a poor  initial  model.  After  a number  of  interations,  PEBLS  had 

reduced  the  RMS  deviation  to  .023,  as  shown  in  the  lower  half  of  Table  XIII. 

The  final  values  of  the  parameters  P.  . , d.  . , E . and  g.  (where  the 

lin'  lin'  plast,  1 

latter  now  represents  the  load-mass  factor,  X)  are  given  in  Table  XIV. 

Again,  it  is  felt  that  these  parameter  values  are  "physically  reasonable" 
and  roughly  in  agreement  with  intuition. 

Biggs'  text,  P.209,  gives  load-mass  factors  in  the  range 

.66  < X < .78 

for  simply-supported  one-way  slabs.  These  values  are  in  contrast  with  the 
factor  1.18  estimated  by  PEBLS  and  given  in  Table  XIV.  Also  note  that  the 
estimated  pressure-vs-def lection  curve  does  not  fit  the  static  data  very 
well  (Figure  5-23) . 

In  retrospect  it  is  thought  that  the  initial  values  (?lin  = Ho, 

d, . = .34)  were  "poor  estimates"  and  that  they  contributed  to  the  relatively 

lin 

poor  lerformance  of  the  estimation  program  in  this  instance.  For  example, 

referring  to  Figure  5-23,  one  sees  that  a better  "initial  estimate"  might 

have  been  (P,  . = 110,  d. . = 1.0)  for  example.  It  is  speculated  that  a 
1 in  lin 

better  initial  estimate  would  have  produced  more  realistic  results  for  the 
load-mass  factor  and  for  the  pressure  vs-def lection  curve. 

* The  reader  will  recall  that  the  "multi-parameter"  model  (Section  5.4) 

selected  a such  that  the  initial  pressures  were  (in  effect)  "cut  off" 

o 

(i.e.,  the  pressure  pulse  was  reshaped  in  time). 
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5.9  DISCUSSION  OF  RESULTS 

These  results  (especially  those  in  Section  5.6,  5.7,  and  5.8)  led  to 
the  conclusion  that  realistic  parameters  can  be  found  that  allow  single  de- 
gree-of-f reedom  models  to  accurately  reproduce  dynamic  response  data  for 
actual  slabs.  It  also  appears  that  (in  many  cases)  these  model  parameters 
are  close  to  those  which  would  be  predicted  a priori. 

A major  result  of  this  effort  is  the  recommendation  that 

• "Tests  on  any  non-linear  structure  (or  element  such  as  rein- 
forced concrete  slab)  should  be  conducted  at  several  levels 
of  force  excitation" . 

When  put  in  the  above  terms,  such  a recommendation  seems  almost  tri- 
vial. Any  student  of  non-linear  vibrations  or  non-linear  circuit  theory 
knows  that  the  dynamic  response  (i.e.,  output)  depends  in  a non-linear  fa- 
shion upon  the  excitation  (i.e.,  upon  the  input).  This  non-linear  depen- 
dence is  the  very  basis  for  labeling  the  problem  as  "non-linear"  in  the 
first  place.  Nevertheless,  such  a non-linear  dependence  between  the  dynamic 
response  and  the  applied  over-pressure  p(t)  was  not  examined  in  two  of  the 
three  main  references  used  herein  as  sources  of  dynamic  test  data.  This 
situation  can  perhaps  best  be  illustrated  by  an  analogy  with  earthquake 
engineering,  as  follows. 

Suppose  one  has  a single-story  steel-framed  structure,  where  the 
behavior  of  the  steel  columns  is  thought  to  be  represented  by  a bi-linear 
(elasto-plastic)  spring.  The  structural  engineer  wishes  to  identify  the 
non-linear  characteristic  (i.e.,  force-deflection  curve)  of  his  structure, 
so  he  can  calculate  if  the  structure  will  fail  during  a major  earthquake. 

If  the  structure  is  excited  only  at  small  amplitudes  (e.g.,  by,  say,  magni- 
tude 3 or  4 earthquakes)  it  will  never  (or  very  seldom)  be  driven  into  the 
non-linear  range,  and  there  will  be  no  way  humanly  possible  for  the  engineer 
to  identify  (with  any  reasonable  degree  of  certainty)  the  non-linear  charac- 
'■•-■■''sties  of  his  structure. 

A similar,  but  less-obvious  condition  exists  when  the  structure  exper- 
iences only  one  non-linear  event,  say  a magnitude  of  6.5  earthquake.  Now 
the  engineer  finds  that  his  structure  has  been  driven  into  the  non-linear 
range,  and  he  can  calculate  (i.e.,  identify)  the  yield  stress,  say,  of  his 
bi-linear  elasto-plastic  force-deflection  curve.  However,  he  still  cannot 
predict  (with  very  much  certainty)  whether  or  not  his  structure  will  collapse 
if  it  is  hit  by  a major  earthquake  (e.g.,  magnitude  7.8,  say).  He  needs 
test  data  at  several  levels  of  excitation  for  his  non-linear  structure. 
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To  return  to  the  problem  at  hand,  namely  the  response  of  R/C  slabs. 
Watt's  tests  (reference  3)  are  analogous  to  the  magnitude  6.5  earthquake 
just  considered.  The  slabs  were  driven  into  the  non-linear  range  (with  a 
maximum  dynamic  deflection  of  about  .4  inches)  but  the  static  load-deflec- 
tion curve  showed  the  slab  still  yielding  at  a deflection  of  .8  inches. 
Without  further  testing,  one  has  little  or  no  way  of  knowing  whether  these 
deep  slabs  will  withstand  a dynamic  overpressure  of  2000  psi  or  20,000  psi 
before  failing  catastrophically.  They  have  not  been  tested  to  a magnitude 
7.8  (major  earthquake  event). 

Conversely,  all  of  Brown  and  Black's  tests  (reference  5)  drove  the 
test  slabs  to  major  failure.  They  measured  a static  failure  deflection  of 
about  4 inches,  yet  each  of  their  test  slabs  was  driven  to  dynamic  deflec- 
tion on  the  order  of  7 inches  or  more.  They  tested  to  only  major  earth- 
quakes (using  our  analogy)  and  experienced  no  moderate  or  low-level  loading 
Clearly,  if  the  structure  is  always  driven  to  collapse,  one  obtains  little 
information  about  its  behavior  throughout  its  dynamic  range. 

This  state  of  affairs  leads  to  the  following  major  recommendation: 

• "Future  dynamic  testing  (of  slabs,  box-like  bunkers,  or  other 
non-linear  structures)  should  be  conducted  at  three  fairly  distinct 
levels  of  excitation.  For  want  of  a more  precise  definition,  the 
three  levels  recommended  are: 

(i)  moderately  non-linear  (e.g.,  2 or  3 elastic  deflections) 

(ii)  strongly  non-linear  (e.g.,  7 or  8 elastic  deflections) 

(iii)  to  complete  failure  (e.g.,  total  collapse) 

Such  a test  sequence  may  (possibly)  triple  the  cost  of  a typical  test  pro- 
gram. However,  an  incomplete  or  inconclusive  test  program  may  prove  mis- 
leading and  might  actually  end  up  costing  more  in  the  long  run.  Keenan 
(reference  6)  recognized  the  need  to  test  at  various  dynamic  pressure 
levels,  which  he  accomplished  with  the  adverse  side-effect  of  having  to  re- 
test partially  damaged  slabs.  Presumably,  Keenan  was  under  a budget  con- 
straint which  limited  the  number  of  virgin  test  specimens  at  his  disposal. 

A difficulty  with  Keenan's  work  (reference  6)  was  that  his  static 
tests  were  not  taken  to  complete  failure  of  his  test  slab.  His  dynamic 
tests  covered  the  ranges  of  linear,  moderately  non-linear,  and  strongly 
non-linear  tests,  but  the  failure  load  (i.e.,  total  collapse)  of  the  slabs 
was  not  determined.  This  situation  led  to  the  following  recommendation: 

• "Static  tests  of  non-linear  structural  elements  (or  complete 
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structures)  should  be  continued  to  the  point  of  total  collapse 
if  at  all  possible." 

The  preceeding  remarks  and  recommendations  are  among  the  major  conclu- 
sions which  kept  recurring  to  the  writers  of  this  report.  They  have  been 
emphasized  herein  with  the  hope  that  they  can  be  utilized  in  planning  and 
conducting  future  tests.  It  is  recognized,  however,  that  in  any  test  pro- 
gram compromises  must  be  made,  and  test  planners  may  choose  not  to  accept 
all  of  our  recommendations. 

The  estimation  program  was  successful  in  fitting  Watt's  deep  slab 
data  with  both  a "multi-parameter"  model  and  a simplified  four-parameter 
model  where  the  values  of  the  parameters  obtained  were  "physically  reaso- 
nable" (i.e.,  they  agreed  with  intuition).  For  Brown  and  Black's  conven- 
tional slab  data,  it  was  again  possible  to  fit  the  dynamic  response  results, 
but  the  result  was  a "near-linear"  model  which  was  much  less  satisfying 
intuitively.  Similarly,  it  was  possible  to  fit  Keenan's  laced  slab  data 
(with  the  multi-parameter  model)  but  the  estimation  scheme  greatly  de- 
emphasized  the  initial  pressures  (of  the  pressure-time  history) . 

Successful  attempts  were  also  made  to  fit  the  various  test  data  with 
simple  three  and  four-parameter  models.  These  results  were  significant,  in 
that  the  models  thus  estimated  were  physically  realistic,  (and  still  quite 
simple)  yet  they  gave  a good  "fit"  to  the  dynamic  data  (i.e.,  low  RMS 
error  values) . Thus,  the  authors  were  led  ultimately  to  the  major  conclu- 
sion that,  "Simple  models  can  be  made  to  work",  as  discussed  in  detail  in 
Section  5.0. 

Regarding  recommendations  for  future  work  in  this  area  of  dynamic 
modeling,  it  may  be  said  that  greater  emphasis  should  be  placed  on  finding 
parameters  which  agree  with  "physical  intuition"  or  elementary  theory, 
rather  than  on  making  a general  model  which  can  encompass  all  the  test 
data.  When  the  study  was  begun,  it  was  not  known  whether  a simple  (3  para- 
meter, say)  model  would  suffice  or  not.  Watt's  report  showed  that  the  slab 
went  through  three  distinct  phases  of  deformation,  and  orginally  it  was 
decided  to  model  deformation-dependent  properties  for  the  mass,  stiffness, 
and  applied  force.  The  result  was  that  (by  using  all  the  parameters)  one 
could  "fit"  the  response  data  rather  easily.  But,  the  problem  (which  did 
not  become  evident  until  later)  was  that  the  parameter  values  obtained  were 
not  readily  acceptable  from  a physical  standpoint.  The  mathematics  of 
system  identification  was  originally  over-emphasized  with  respect  to  the 
physics  of  the  slab-deformation  problem. 


It  is  also  recommended  that  some  "sensitivity  studies"  be  conducted, 
namely  to  see  how  much  that  maximum  displacement  (Y  max)  veries  with  changes 
in  the  parameters  (A,  Plin»  diin'  etc-)  For  example,  just  how  significant  is 
the  fact  that  the  load-mass  factor  was  estimated  at  .82  (vs.  a theoretical 
value  ranging  from  .51  to  .67)  for  Brown  and  Black's  slabs.  Would  a value 
of  A = .7  have  changed  the  results  much?  The  answer  to  these  questions  are 
left  to  a future  study. 


6. 


CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


The  primary  conclusion  of  this  report  is  that 

. Realistic  parameters  can  be  found  (using  Parameter  Estimation) 
that  allow  single  degree-of-freedom  models  to  accurately 
reproduce  dynamic  response  data  for  actual  slabs. 

A related  secondaty  conclusion  is  that 

Often  (but  not  always)  these  model  parameters  can  be  estimated 
fairly  accurately.  The  estimates  can  be  provided  by  static  test 
(if  available) , and  by  yield-line  theory  if  test  data  is  missing. 

These  conclusions  are  substantiated  in  the  body  of  the  report. 

Major  recommendations  of  this  report  are  that 

• "Tests  on  any  non-linear  structure  (or  element  such  as  reinforced 
concrete  slab)  should  be  conducted  at  several  levels  of  force 
^vci tation" . 

• "Future  dynamic  testing  (of  slabs,  box-like  bunkers,  or  other 
non-linear  structures)  should  be  conducted  at  three  fairly 
distinct  levels  of  excitation.  For  want  of  a more  precise 
definition,  the  three  test  levels  recommended  are: 

(i)  moderately  non-linear,  (e.g.,  2 or  3 elastic  deflections) 

(ii)  strongly  non-linear,  (e.g.,  7 or  8 elastic  deflections) 

(iii)  to  complete  failure  (e.g.,  total  collapse) 

• "Static  tests  of  non-linear  structural  elements  (or  complete 
structures)  should  be  continued  to  the  point  of  total  collapse 
if  at  all  possible." 

The  rationale  behind  these  recommendations  is  included  in  the  main 
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Least  Squares 


In  Eq.  2,  the  vector  x^  represents  the  independent  variables 
used  to  generate  the  calculated  response  and  for  a ground 
shock  calculation  are  represented  by  location  and  time.  Of 
crucial  importance  to  the  least  squares  method  is  that  the 


equation  chosen  for  u^  is  the  best  representation  of  the  experi- 
ment being  simulated. 

By  the  Gauss  or  least  squared  method,  we  seek  to  minimize  the 


sum  of  the  squares  residuals  given  by 


where  n is  the  total  number  of  experimental  measurements. 

Thus,  Eq.  (3)  is  the  objective  function  first  evolved  by  Gauss. 
Note  that  this  objective  function  does  not  account  for  pos- 
sible errors  in  the  independent  variable.  Experience  has  shown 
that  errors  in  the  dependent  variable  are  usually  two  orders 
of  magnitude  higher  than  those  in  the  independent  variable 


thus  allowing  the  latter's  errors  to  be  neglected.  Nonethe- 
less, fo?  those  special  cases  where  the  two  types  of  errors 
are  of  ecual  magnitude.  Eg.  (3)  can  be  suitably  modified  (see 
Ref . 1)  . 


If  the  measurements  x.  and  u are  repeated  n times  and  the 

1 °i 

average  values  computed  these  average  values  should  approach 
the  true  values  as  n approaches  infinity.  This  will  occur 
if  systematic  errors  do  not  affect  the  measurements.  Such 
errors  consistently  cause  the  measured  response  to  be  either 
larger  or  smaller  than  the  true  values.  The  errors  can  arise 
from  three  sources: 

1.  A mistake  in  the  measuring  technique 

2.  A fault  in  a measuring  instrument 

3.  Failure  to  account  for  a factor  affecting 
the  experiment. 

Tnus  for  example,  if  temperature  affects  the  response  of  an 
accelerometer,  such  a devise  located  in  close  proximity  to 
the  heat  of  a blast  might  contain  a systematic  error  in  its 
measurement. 

The  simplest  case  of  least  squares  is  that  of  the  polynomial. 
For  a straight  line  this  becomes 


1 1 

- Z <V  • rl  ‘ r2xi> 


where  for  the  sake  of  simplicity  we  are  considering  the  case 
of  a single  independent  variable,  x^.  We  see  than  that  the 
objective  function,  F,  is  a function  of  r^  and  r^.  Of  the 


many  values  that  these 
set  that  minimizes  F. 


parameters  can  assume,  we  seek  the 
At  this  minimum,  we  require  tha\ 


(5) 


(6) 


which  yields  two  equations,  (It  should  be  noted  that  no 
maximum  exists  for  F since  for  any  maximum  one  might  postu- 
late, an  even  greater  maximum  can  be  found.)  Applying  the 
criteria  of  Eqs.  (5)  and  (6)  to  Eq.  (4)  yields  the  two  linear 
equations 


rln 


+ r. 


r 


1 


x . u 
1 o . 
1 


(7) 


(8) 


which  can  be  readily  solved  for  and 

The  above  derivation  can  be  extended  to  any  polynomial  of  order 
p-1  defined  by 


k-1 


x . 
r 


(9) 


If  we  cast  the  resulting  set  of  equations  (7)  and  (8)  in  matrix 
form  obtaining 

Cr  = V (10) 
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it  can  readily  be  shown  that  the  el  .-ments  of  C and  V are 
given  by: 


and 


V . 
D 


n 


(ID 


( 12) 


An  interesting  phenomenon  arises  when  the  degree  of  the 
polynomial  becomes  too  high.  If  we  normalize  our  independent 
variable  so  that  it  ranges  from  0 to  1 (which  is  always  pos- 
sible) and  employ  a sample  of  uniformly  distributed  data, 
then  the  elements  of  C are  given  by 


'jk 


_ f1  j+k-2j  n 

* X x dx  = 3+k-r 


(13) 


The  terms  represented  by  are  elements  of  the  so-called 

Hilbert  matrix.  The  determinant  of  this  matrix  is  given  by  ( 


n = [H213J (P-1)  ! ] 3 

P p* (P+1)  1 . . . (2P-1)  1 


(14) 


which  approaches  zero  very  rapidly  as  the  matrix  grows  as 
shown  in  Table  1. 





p 

DP 

1 

2 

e.sx’D'2 

3 

4.6X10'4 

4 

1 .7X1 Cf 7 

5 

-12 

3.7X10 

6 

5 .4X1 0"1 8 

7 

4.8X1 0-25 

8 

2 . 7X1 0- 33 

9 

9. 7X1 0"43 

Table  1.  Determinants  of  Hilbert-like  Matrices  as  a 
Function  of  Matrix  Size 

When  the  determinant  of  the  matrix  associated  with  a linear 
set  of  equations  approaches  zero,  the  system  is  said  to  be 
ill-conditioned . The  set  of  parameters,  r,  generated  from 
such  a system  are  adversely  affected  by  round-off,  i.e.  it 
has  a large  error  associated  with  it.  In  the  extreme,  the 
determinant  will  be  so  close  to  zero  so  as  to  render  the 
matrix  to  be  computationally  singular.  Experience  has  shown 
that  this  occurs  when  the  number  of  parameters  is  about  6 
(fifth-order  polynomial).  The  method  most  often  used  to 
circumvent  this  problem  is  the  implementation  of  orthogonal 
polynomials.  However,  since  the  desired  effect  is  somewhat 
artificial,  it  is  preferable  to  employ  a lower  order  poly- 
nomial (fourth-order  or  lower). 


LJ 


The  least  squares  method  is  not  confined  to  polynomials.  The 
method  can  easily  be  extended  to  other  functions  by  allowing 
u^  to  assume  a form  that  is  linear  with  respect  to  the  para- 
meters : 

ui  = rl9i^xi)  + r2g2(xi)  + + r pg p(xi> 


i. 

~ 2.  r|cW 


k=l 


(15) 


By  inserting  Eq.  (15)  into  Eq.  (13)  and  requiring  that 
9F 


3r, 


= 0 k = 1 , 2 . . . p 


(16; 


we  again  obtain  a linear  set  of  equations  of  the  form 


Cr  = V 


(10) 


where  now  the  elements  of  C and  V are: 


n 

- -i  = / g . (x  . ) g,  (x  . ) 

jk  i'  i 


(17) 


and 


n 


v* = ^ V*i)u°i 


(18) 


It  is  seen  that  the  polynomial  of  Eq.  (9)  is  a special  case 
of  Eq.  (15)  where 


/ , k-1 

gk(xi}  = xi 


(19) 


Since  most  functions  are  expandable  into  polynomials,  the 
earlier  comments  regarding  the  Hilbert-like  matrix  generated 
by  the  use  of  high-order  polynomials  can  be  extended  to  most 
functions  defined  by  Eq.  (15) . 
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A cursory  examination  of  C_.k  in  Eq.  (17)  will  reveal  that 
C is  a symmetric  matrix.  This  property  which  also  app  ies 
to  the  matrix  generated  through  the  polynomial  (see  Eq  (11)) 
permits  a substantial  savings  in  computational  time  si.;ce  only 
the  upper  or  lower  triangle  of  C need  be  computed  (al  >ng  with 
its  diagonal  elements) . 


For  ease  of  notation,  define  the  matrix  T by 

Ti j = 9 j (*i ) (20) 

where  T is  n * p.  Then  in  matrix  form,  the  parameter  estimate, 
r,  in  Eq.  (10)  is  given  by 

T - 1 T 

r = (T  T)  T u • (21) 

o 


Therefore,  the  matrix  given  by 

G = (TTT)_1  TT  (22) 

is  that  vehicle  which  yields  the  estimate,  r,  as  a linear 
transformation  of  the  experimental  response  (vector  of 
observations) . 

It  is  obvious  that  the  least  squares  technique  cannot  be 
employed  when  the  number  of  parameters  exceeds  the  number 
of  measurements,  i.e.  p > n.  When  this  occurs,  the  system  of 
Eq.  (10)  is  said  to  be  underdetermined . When  p is  equal  to  n 
each  of  the  residuals  in  Eq.  (3)  will  be  zero  resulting  in  a 
minimum  value  of  zero  for  F.  For  this  case  the  system  is  said 
to  possess  zero  degrees  of  freedom  which  is  the  difference 
between  p and  n. 
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Weighted  Least  Squares 
The  formulation  of  Eq.  (3) 


n 


has  two  implicit  assumptions:  firstly,  that  all  the  experi- 
mental points  are  measured  with  the  same  degree  of  precision, 
and  secondly,  that  all  the  residuals  have  the  same  dimensions. 
Therefore,  if  some  points  were  measured  with  high  precision, 
Eq.  (3)  would  treat  these  points  with  equal  weight  as  the 
other  observations.  As  an  example  of  the  second  assumption, 
if  some  of  the  responses  were  in  units  of  seconds  and  the 
rest  in  units  of  meters/second  (as  might  happen  in  a velocity 
waveform  response  to  a ground  shock)  the  formulation  of  Eq. 

(3)  would  be  untenable. 


The  above  two  problems  are  solved  by  adding  statistical  weight- 
ing to  the  least  squares  formulation.  Each  squared  residual, 

2 

U^,  is  multiplied  by  a weight, w^ , which  is  equal  to  the  inverse 
of  the  square  of  the  uncertainty  or  standard  deviation  or: 


*“■77 

l 


(23) 


c e 


ii 


t h 

where  S£r.,  is  the  variance  of  the  i observation.  The  re- 

ii 

suiting  weighted  least  squares  formulation  is  given  by 


- I ...  o2 

11  1 


mm 


(24) 


i=l 
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Minimizing  this  scalar  with  respect  to  each  parameter  again 
yields  a linear  system  of  equations: 

Cr  = V (10) 


where 


C 


jk 


wiigj(xi)gk(xi) 


(25  ) 


and 


V,. 


11 

= I V.. 

11 


i=l 


VVV 

1 


(26  ) 


It  should  be  noted  the  C matrix  defined  through  Eq.  (25)  is 
again  symmetric. 

In  matrix  form,  the  parameter  estimate,  r,  in  Eq.  (10)  is 
now  given  by 

T -IT 

r = ( T wT ) T wuQ  (27) 


where  the  matrix,  T,  is  again  given  by: 


T.  . 

il 


g.(xi) 


and 


(20) 


(23) 


where  both  w and  are  diagonal  matrices.  Thus,  the  linear 

transformation  from  u to  r is  achieved  by 

o 

r = Gu  (28) 

o 


A-ll 


where 


T - 1 T 
G = (T  wT)  T w 


Now  the  covariance  of  the  estimate  is  given  by 


* * * t 

£rr  E E(ZrrZrr» 


Zrr  is  the  random  error  in  r,  and  E denotes  expectation.  By 
the  transformation  of  Eq.  (28): 

Srr  ’ E[GZ„IGZtE)T]  = GEUeczJ)GT  (31) 

T 

where  E(Z£;rZ_c)  is  the  expectation  of  the  observation  error, 
u - u , and  related  to  S by 

See  = E(ZccZee}  (32) 

where  Z is  the  random  error  in  u . 

sc  o 

Returning  to  Eq.  (31) 

if  T T T 

S = GE ( Z Z )G  = G S G1 

rr  cc  cc  ec 


( TTwT) ” 1 TTww-1 [ (TTwT) _1  TTw]T 


( TTwT ) ~ 1 TT[(TTwT)"1  TTw]T 


(TTwT)  "-1  [ (TTwT)  "1  (TTwT)  ]T 


( TTwT)  ” 1 = (TTS"3T)_1=  c-1 


At  this  point  we  should  now  pose  the  question,  "How  do  we 
know  that  the  estimate  of  the  experimental  error  of  the  response 
is  correct?"  This  estimate  is  based  on  the  judgement  of  the 


experimentalist  and  when  the  measurements  are  conducted  by  a 
team  it  is  unusual  for  a group  to  agree  on  an  error  estimate. 

To  answer  the  above  question,  we  resort  to  Eq.  (32) 


Having  chosen  a mathematical  model  for  our  experiment,  we 
believe  this  model  to  be  correct.  Thus,  if  there  were  no 
experimental  errors,  and  the  correct  model  was  postulated, 
our  experimental  points  would  coincide  with  those  of  the 
model,  i.e.,  all  the  residuals  would  be  zero.  Therefore, 
the  residual  Z^r  appearing  in  Eq.  (32)  can  be  attributed 
to  experimental  and  mathematical  error  and  can  then  be  used 
as  a check  on  the  input  covariance.  If  Eq.  (32)  does  not 
yield  a covariance  matrix  similar  to  the  input  covariance 
matrix,  then  S should  be  adjusted  and  the  least  squares 
process  repeated.  When  Eq.  (32)  is  finally  satisfied,  the 
error  estimates  will  be  substantially  correct.  Unless  the 
above  steps  are  carried  out,  Eq.  (33)  , v;hich  specifies  the 
revised  parameter  covariance  will  be  incorrect. 

When  the  above  errors  are  experimental,  we  have  implicitly 
assumed  that  the  mathematical  model  is  correct.  When  the 
model  is  incorrect,  we  have  failed  to  account  for  factors 
affecting  the  experiment  and  possibly  introduced  a systematic 
error.  Under  this  condition,  the  residual  ZC{.  in  Eq.  (32) 
will  not  be  completely  generated  by  experimental  error.  How- 
ever, it  is  reasonable  to  assume  that  if  several  models  are 
postulated  and  the  above  steps  executed,  such  that  Eq.  (32) 
is  obeyed,  then  the  model  with  the  best  fit  would  be  closest 
to  the  true  mathematical  model. 
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The  technique  of  adjusting  the  covariance  matrix,  S , 

c c 

should  not  be  interpreted  as  casting  doubt  on  the  credibility 
of  the  experimentalist  who  first  estimates  the  experimental 
errors.  If  the  experimentalist  is  regarded  to  be  competent, 
then  a substantial  disparity  between  the  original  and 
adjusted  covariances  would  most  likely  indicate  that  the 
mathematical  model  is  incorrect  and  a reappraisal  is  suggested. 


Nonlinear  Weighted  Least  Squares 


T.ne  case  of  nonlinear  least  squares  occurs  when  the  function 
appearing  in  the  objective  function 


F = 


n 

^ * 
1=1 


1 1 

2 y , 

iiL’i  = Zwn(u 


i=l 


' V 


(34; 


is  a nonlinear  combination  of  the  parameter  vector  r: 


U.  = u (xif  r , r ...  r ) (35) 

Note  that  for  the  sake  of  clarity,  we  have  assumed  only  one 
independent  variable  in  Eq.  (34).  Nonetheless,  the  ensuing 
derivation  can  readily  be  extended  to  multiple  independent 
variables  by  substituting  a vector  for  x^. 

The  reader  can  readily  determine  that  minimizing  Eq.  (34) 
using  Eq.  (35)  leads  to  a set  of  nonlinear  equations.  Rather 
than  treat  each  case  individually,  it  is  preferable  to  start 
with  Eq.  (35)  and  develop  a general  technique. 


Minimizing  the  objective  function  of  Eq.  (34)  with  respect 
to  each  parameter,  r^ , yields 


n 

V 


9u 


W . . (U  - U . ) -r— 

i<L  11  °i  1 9rk 


i = 0 


k = 1,  2,  3 . . . p (36) 


In  order  to  obtain  a linearized  formulation,  we  expand  u into 
a truncated  Taylor  series  evaluated  at  an  estimated  value  of 
r denoted  by  the  index,  e: 


p 

jI1 


(15) 


then 


gk(xi> 


(41) 


Therefore  the  T matrix  defined  by  Eq.  (38)  is  the  same  as  the 
T matrix  specified  by  Eq.  (20). 


« 
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Returning  tc  che  n 'nlinear  problem,  we  take  the  derivative 
of  in  Eq . (40)  ith  respect  to  r,  : 


c U . 

^k  " Tik 


(4.2) 


Inserting  Eqs.  (40)  and  (42)  into  Eq . (36)  yields 

\ i ■ \T  . , 

P ip/  ik 


4WiiK  ' % - Vil  - Vi2 


® k - 1 , 2 , 3 . . . p 

The  expression  of  Eq.  (43)  gives  the  linear  set  of 
of  the  form 

CR  = V 

where  the  elements  of  C and  V are 


Cjk 


4 i 

■ S w 


i=l 


• T.  .T., 

i 13  ik 


(43) 

equation: 

(44) 


(45) 


and 


n 

* 


(46) 


Note  that  C is  symmetric. 

When  the  system  of  Eq.  (44)  is  solved  for  R,  the  new  estimate 
ot  r is  given  by 


r = r + R 
e 


(47) 
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This  value  of  r is  only  approximate , but  if  the  sequence  is 
converging,  it  is  a better  estimate  than  the  previous  set 
in  r . Therefore,  if  we  reolace  the  values  in  r.  by  those 
from  r and  repeat  the  sequence,  we  will  get  an  even  better 
set  for  r.  In  the  limit  as  r approaches  re , the  higher  order 
terms  in  Eq.  (37)  go  to  zero  faster  than  R,  and  the  formu- 
lation becomes  exact. 

The  choice  of  the  initial  guess  for  re  in  the  above  sequence 
can  be  of  paramount  importance.  Experiments  for  nonlinear 
cases  have  shown  that  some  cases  have  converged  when  the  initial 
and  final  values  differed  by  a factor  of  10^,  while  in  other 
cases,  the  sequence  diverged  when  the  initial  values  dif- 
fered from  the  true  values  by  only  20%.  In  general,  however, 
the  better  the  initial  guess,  the  greater  the  chances  for 
convergence.  When  divergence  does  occur,  other  methods  can  be 
employed.  Often  a seemingly  diverging  problem  can  be  brought 
to  convergence  by  implementing  a search  algorithm  (see  Ref.  3) . 

Eqs . (44)  - (47)  are  quite  general,  and  can  be  applied  to  the 

linear  case.  For  this  purpose,  we  set  r0  equal  to  zero.  The 
reader  can  readily  determine  that  the  linear  case  is  then 
retrieved,  and  is  solved  in  a single  iteration. 


In  matrix  form,  the  parameter  estimate  r in  Eq.  (47)  is 
r = re  + R = re  + C_1V 

= re  + (TTwT)_1TT  w(uo  - ue) 


or 


r = re  + G(uo  - ue) 


(48) 


A-18 


where  G has  been  derived  previously  as 


T -1  ? 

G = (T  wT)  T w 


The  covariances  of  the  estimate  and  the  observation  are 


respectively 


★ * ★ T* 

S = E (Z  Z ) 
rr  rr  rr 


S = E(Z  2r) 

CC  C C C € 


As  has  been  derived  earlier  Z , is  related  to  Zrr  by 


Z = GZ,r 
rr  «-  c 


Here,  however,  we  must  assume  that  the  true  value  of  r and 
the  calculated  value  of  r are  sufficiently  close  so  that 
the  G matrix  is  the  same  for  both.  With  this  assumption, 
we  continue,  and  applying  Eq.  (33)  obtain  the  revised  co- 
variance  matrix. 

* T - 1 - 1 

S = (T  S 1 T)  1 (33) 

rr  cc 

As  has  been  noted  for  the  linear  case,  this  derivation  is 

based  on  the  assumption  that  Eq.  (32)  for  the  observation 

T 

covariance  matrix  is  valid.  If  the  expectation  of  Z__Z.„  is 

not  equal  to  S , then  S should  be  adjusted  and  the  sequence 
^ cc  cc 

of  iterations  repeated  until  we  obtain  a close  approximation 
for  S££  in  E(2c£2^I. 
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The  problem  of  ill-conditioning  associated  with  the  C matrix 
in  linear  least  squares  is  also  found  to  occur  with  the  C 
matrix  of  Eq.  (45)  in  the  nonlinear  mode.  This  is  because 
most  functions  can  be  approximated  by  a power  series.  This 
author  has  found  the  problem  to  be  more  severe  in  the  nonlinear 
case  having  noted  severe  ill-conditioning  commonly  manifested 
when  5 parameters  were  used. 


I 


Minimum  Variance 


Up  to  now  we  have  assumed  that  the  errors  in  the  measu  enents 
are  statistically  independent  so  that  the  covariance  .atrix 
See  is  diagonal.  In  many  experiments  a cross-correlat  on 
exists  such  that  other  elements  of  S are  non-zero.  One 
important  property  of  this  matrix  is  the  fact  that  it  must 
be  symmetric  since  the  matrix  is  defined  by  E(z  Z7  ) . 

EE  EE 

For  the  most  general  case  which  is  nonlinear  the  objective 
function  is 


n n 


or  in  matrix  form 

T 

F = (uQ  - u)  w(uq  “ u)  -►  min  (51) 


where  w is  a weight  matrix  defined  by 


and  since  S is  symmetric 

w = wT.  (53) 

Minimizing  the  objective  function  of  Eq.  (50)  with  respect 
to  each  parameter,  r^,  yields: 


= S J W,^T,U  /u 
i j 


ij*ik(uOj  ue  j ) k “ 1'  2'  3...p  (57) 


The  expression  of  Eq.  (57)  yields  the  linear  set  of  equations 
of  the  form 


CR  = V 


where  the  elements  of  C and  V are  defined  by 

ck£  5 ISwiiTi: 


1 3 


'i j * ikT 


(58) 


(59) 


and 


■ 2S“ijTik(uo  - «e  ) 
i 3 '3  3 / 


(60) 


A-22 


In  matrix  form  the  parameter  estimate,  r,  imbedded  in  Eq . (58) 

is  given  by: 

r = re  + (TTwT)_1  TTv(uo  - u£)  (61) 

As  in  the  case  of  nonlinear  least  squares  the  solution  is 
obtained  by  iteration.  The  matrix  used  for  the  linear 
transformation  is  seen  to  be 


T -1  T -1  T 

G = (T  WT)  T1  w = C T w . 


It  is  interesting  to  note  the  difference  between  least  squares 
and  minimum  variance  as  embodied  in  the  transformation  matrices 
of  Eqs.  (29)  and  (62).  The  weight  matrix  w appearing  in 
Eq.  (29)  is  diagonal  while  its  counterpart  in  Eq.  (62)  is 
not.  However,  as  seen  in  Eq.  (59),  since  the  inverse  of 
C is  symmetric,  C is  also  symmetric  which  is  also  true  for 
the  least  squares  case.  Similarly  the  dangers  implicit  in 
the  ill-conditioning  of  the  C matrix  resulting  from  the  use 
of  too  many  parameters  can  also  be  carried  over  to  minimum 
variance. 

The  covariance  matrix  associated  with  the  parameter  estimate 
is  specified  by  Eq.  (33) 

* T -1  -1  -1 

srr  = (T  S^T)  = C (33) 

It  should  be  noted,  however,  that  S as  defined  in  this 
section  will  now  be  full  and  not  diagonal. 
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Bayesian  Estimation 


As  seen  earlier,  the  technique  of  minimum  variance  cons- 
titutes the  most  general  tool  for  parameter  estimation. 
Theoretically,  thf'  problem  is  well-posed  and  barring  any 
ill-conditioning  in  the  C matrix  and  divergence  for  the  non- 
linear case,  has  a well-defined  solution.  However,  for  a 
substantial  number  of  observations,  the  weight  matrix,  being 
of  size  L x L,  the  storage  capacity  of  the  computer  may  be 
exceeded.  The  problem  can  be  resolved  by  dividing  the 
experimental  data  into  batches  and  processing  them  sequentially. 
However,  to  accomplish  this  we  must  transfer  from  batch  to  batch 
the  knowledge  we  have  gained  concerning  the  parameter  estimate. 
Mathematically,  this  technique  is  known  as  "Bayesian  Estima- 
tion" after  the  mathematician  who  first  expressed  the  concept 
in  terms  of  probability. 


In  Bayesian  estimation  we  are  given  a prior  estimate  of  the 

parameters,  rQ,  along  with  the  associated  covariance  matrix 

S . We  then  seek  to  minimize  the  objective  function 
rr  J 

n n 


= 2*  2^  w . . (u  -ulfu  - u. 
jti^i  i]\  o.  i/\0j  3) 

P P 

+ zL  ^ w . . (r  - r-Wr  - r 


where  w = S 


and  w = S 


and  S is  a symmetric  matrix  as  is  S 


In  comparing  this  objective 


function  with  that  of  Eq.  (50),  we  find  that  the  second  double 


series  term  on  the  right-hand  side  of  Eq.  (63)  accounts  for 
our  knowledge  of  the  Bayesian  pr^or.  Then  logically  the  new 
parameter  estimate  we  will  obtain  via  Eq.  (63)  will  be  a 
compromise  between  our  knowledge  concerning  the  experimental 
data  (the  first  double  series  term  in  Eq.  (63))  and  the  Bayesia 
prior.  This  means  that  were  we  to  solve  for  the  new  parameter 
estimate  using  Eq.  (50),  the  solution  would  be  the  same  as 
that  for  the  extreme  case  for  Eq.  (63)  when  S -+  «. 

Minimizing  Eq.  (72)  with  respect  to  each  parameter  r^  gives 


Since  the  function  u is  nonlinear,  we  expand  it  again  into 
a truncated  Taylor  series  evaluated  at  an  estimated  value,  r&  : 


u*  * ue£  + R1TU  + R2T£2  + 

and  note  that 
9uf 

= T 

9rk  £k 


R T0 
P *P 


(67) 


(68) 
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Substituting  Eqs.  (75),  (76)  and  (77)  into  Eq.  (74)  gives 

^ ^\w../u  - u - R T . . - R 0T  . _ - ...  - RT  W., 

i ~ 17\  °j  ej  1 ]1  2 D 2 P DP/  ik 


+ X X W . . (u  - u,  - R.T.,  - R _T  • _ - ...  - RT.  \t  .. 
i j 13 \ °i  i 1 11  2 12  F lp/  :k 


^ w.  . (r  - r . \ + Z,w../r  - r ■ | = 0 

, M o,  c)  Y - k ^ o,  i} 


Since  Se£  and  Srr  are  symmetric,  w and  w are  also  symmetric. 
Therefore  the  two  double  series  terms  of  Eq.  (70)  are  identical 
as  are  the  two  single  series  terms.  Eq.  (70)  then  simplifies 


? • vn  • 

* ? "k,i  (v  - r.')  = 0 


R_T._  - ...  - R T.  T . . 
2 D2  P DP  ik 


We  then  note  that 


ro.  ‘ r/  = ro.  - re,  - * 


which  upon  insertion  into  Eq.  (71)  and  some  manipulation 
gives 
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Rl(Wkl+  ^ ^WijTikTjl) 

+ r2  |wk2  + ? 4<  wijTikTj2  ) 

+ ...  + R (w,  + ) Zw..T.,T.  ) 

P V kp  ^ j ID  ik  dp  ) 

p n n 

= A w.  . /r  -r  |+A/u  -u  \ 2 w 

& ‘-ri  v v i-i 


■ ■ T . . 
ID  lk 


k - 1 # 2 , ...  p ( "^  u ) 

The  expression  of  Eq.  (73)  gives  the  linear  set  of  equations 
of  the  form 


CR  = V 


where  the  elements  of  C and  V are  defined  by 
n n 


'kZ  ™kZ  ^WijTikTj2 

i=l  d=1  J 


p 11  11 

= I w,  , /r  -r\+^/u  - u \ X 

A *j[od  e.f  j = 1[oj  e.}  i=1 


w . T . , 
ID  ik 


As  seen  from  Eq.  (75) , the  C matrix  is  symmetric  which  again 
eases  much  of  the  computation. 

In  matrix  form  the  parameter  estimate  is  obtained  by  iteratively 
solving  the  relation 


r = r + 
e 


w + TTwT  ) ” 1 w (r  - r + TTw  (u  - u ) i (77) 
\ o e V o e M 
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The  iteration  starts  with  r replacing  re,  calculating  T as 
a function  :>f  the  independent  variables  and  r , and  then 
generating  : by  Eq.  (77).  When  r converges  to  re , the  itera- 
tions are  1 ^rminated. 


An  alternate  expression  for  r will  now  be  determined.  If  the  term 
(ro  ~ re  ) <7e> 


T 

T WT! 


is  added  and  subtracted  to  the  terms  appearing  in  the  square 
brackets  on  the  right-hand  side  of  Eq.  (77) , we  obtain 


r = re  + 


(w  + TTwt)  -1  w(ro  - re) 

tTwT  (ro  ' re)  - TTwT(ro  - rj 
tTw(uo  - ue) 


re  + (w  + T^wT j~^w  + TTwT)(rQ  - rg) 

|w  + TTwT  j ~ ^ -TTwT  {io  - re)  + TTw(uo  - ug) 
re  + (ro  - re) 

(w  + TTwt)_1  TTw(uo  - ue)  - TTwT (rQ  - rj 


(; 


T \ - 1 T 

= rQ  + | w + T wT ) T w 


= r +Gu  -u^-Tr  - r 
o o e \ o e 


u - u - T ( 

o e \ 

lrc  - re  ) 


r - r 
o e 


)] 


(79: 
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G 


(80) 


■ 


where 


E (w  + TTv  r ) -1TTw 


The  revised  covarieice  matrix  of  the  parameter  is  found  by 
noting  that  the  err  ar  of  the  new  estimate  is  related  to  the 
error  in  the  observation,  zec'  and  the  prior  parameter  error 
Zrr  by  Eq.  (77) 

T 

+ T wT  , 


:*r  = (w  + tTwt j~^w  Zrr 
+ (w  + TTwT ) ” 1 TTw  Z 


ec 


(81) 


The  revised  covariance  matrix  is  then 

* I * * t \ 

S = E Z Z 
rr  1 rr‘rr  / 

= (w  + T^t)*1^  (zrrZ 
+ (w  + TTwT | ~ 1 TTwE | 


T 

rri 


) (w  + TTwt|_1w 


T 

z z 

ee  ee, 


| |w  + TTwt|_1TTw 


(82) 


where  the  cross-correlation  between  Zrr  and  Zec  is  assumed 
to  be  zero.  Noting  that 

•T.)  -c-1 


E zrrzrr, 


(83) 


and 


e(z,ezL)  - 


(84) 


gives  upon  substitution  into  Eq.  (82) 
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rr 


(w  + TTwt)”'lw_1  [ (w  + TTwt)-1w 

T 

(w  + TTwt)'1TTwv~1  t (w  + TTwt)_1TTw] 

T 

(w  + TTwt)_1[(w  + TTwt)-1w] 

rTwT)_1  (w  + ttwt)-1ttwt] 
^wt)"1  (w  + TTwt)_1(w  + TTwt) j 
'l'1  = C'1 


(w  + 
(w  + 
(w  + 


T wT 


(83) 


As  in  case  of  weighted  least  squares  Eq.  (82)  assumes  that 
the  true  value  of  r and  the  calculated  value  of  r are 
sufficiently  close  so  that  the  T matrix  is  the  same  for 
both.  If  this  is  not  true,  a considerable  error  may  be 
introduced . 


To  summarize  the  technique  of  Bayesian  estimation,  we  itera- 
tively solve  the  sequence 


T = T 


(re ) 


(84) 


S = 


\w  + 


T wT 


1 ■ 


(85) 


G = ST 


rw  = G (re) 


= r + G u - u - T ( 
o [ o e ' 


r_  - r. 


(86) 

(87) 


and  repeat  the  next  sequence  with  r replacing  rg.  We  con- 
tinue to  iterate  until  the  difference  between  r and  r 


A-30 


becomes  sufficiently  small  and  the  sequence  has  therefore 
converged.  The  final  value  of  S is  S*r,  the  revised  para- 
meter covariance  matrix.  It  should  be  noted  that  when  the 
response  is  linear  with  respect  to  the  parameters,  only 

one  iteration  is  reauired.  This  means  that  since  r is 

e 

always  set  equal  to  r^  for  the  first  iteration,  the  term 

(r  - r ) in  Eq.  (87)  is  zero  for  the  linear  case.  Alter- 
o e 

natively,  as  will  be  noted  later,  one  can  choose  to  set  both 
re  and  ue  equal  to  zero  for  the  linear  case. 

When  sequential  batch  processing  is  used,  the  Bayesian  tech- 
nique should  yield  the  same  solution  as  that  which  would 
be  obtained  by  processing  all  batches  together  without  a 
Bayesian  prior.  To  test  this,  let  us  confine  ourselves  to 
two  batches  denoted  by  indices  1 and  2 along  with  a linear 
estimator.  The  processing  of  the  first  batch  uses  no  prior 
estimate  and  gives  a parameter  estimate  (Eq.  (77)): 


= (tjVi)-1  T1Tw1uc 


and  a revised  covariance 


Srrj  ■ KVl)'1  * V1  (89> 

The  estimate  and  the  covariance  of  Eqs.  (88)  and  (89)  are 
input  as  information  about  the  prior  for  the  second  batch. 
This  gives  (Eq.  (77)) 

r2  ‘ ("2  + T2T“2T2)"1(“2rl  + T2Tw2uoJ 

* (TlTVl  + T2T“2T2)'1(T1Tw1uo1  + T2T“2uo2) 
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Ea . (90)  gives  the  estimate  resulting  from  a linear 

accumulation  of  all  the  experimental  data,  hence,  proving 
the  affect  of  the  Bayesian  prior.  For  K batc'ies  of  data 
this  gives 


T 

T . w . T . 
l 11 


V1  y t 

I ? T.  w.u 
. I 1 1 o . 

/ 1=1  1 


The  proof  can  readily  be  extended  to  the  nonlinear  case. 

Here,  however,  the  final  value  of  T.  for  each  batch  will 
not  necessarily  be  the  same  thereby  yielding  a different 
estimate  for  the  final  batch.  If  the  system  defining  u is 
highly  nonlinear  such  that  large  variations  in  are  en- 
countered from  batch  to  batch,  the  error  could  be  substantial. 

The  procedure  outlined  in  this  section  is  correct  provided 
that  the  covariance  matrices  of  the  observations  and  the  prior 
are  correct.  As  noted  in  the  earlier  sections,  the  new 
parameter  estimate  has  little  merit  if  wrong  covariances  are 
employed . 

Eqs.  ( 84 ) through  (87)  constitute  the  most  general  case  of 
parameter  estimation,  and  the  reader  will  find  that  the 
equations  of  the  previous  sections  are  special  cases  of 
this  system.  However,  the  Bayesian  aspect  which  has  been 
added  in  this  section  is  particularly  appealing  in  the  non- 
linear case  for  three  reasons: 

1.  The  prior  parameter  estimate  constitutes  a good  guess 

of  the  parameter  estimate  and  hence  increase  the  chances 
for  convergence. 
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2.  The  Bayesian  formulation  drives  the  solution  closer  to 
the  prior  than  does  minimum  variance  (the  Bayesian  :ech- 
nique  with  no  prior).  This  also  leads  to  faster  cc  - 
vergence . 

3.  As  seen  from  Eq . (69)  the  dependency  of  the  Bayesian 

formulation  on  the  prior  is  linear.  This  means  that  as 
the  dependency  on  the  prior  is  increased  (via  S^)  the 
problem  becomes  more  linear  and  convergence  is  enhanced. 

In  addition  to  the  above  observations,  the  addition  of  a 
Bayesian  prior  to  both  linear  and  nonlinear  problems  renders 
the  C matrix  of  the  resulting  set  of  equations  more  diagonally 
dominant  and  hence  more  non-singular.  This  is  why  experiments 
have  shown  that  when  a prior  estimate  is  used,  the  number  of 
parameters  may  be  increased  substantially  without  generating 
a Hilbert-like  matrix. 

The  use  of  a prior  estimate  is  not  restricted  to  sequential 
batch  processing.  Often,  for  example,  in  ground  shock  calcu- 
lations, the  analyst  may  want  to  process  a single  batch  of 
experimental  data  and  incorporate  what  he  believes  to  be  a 
reasonable  estimate  of  the  parameters,  along  with  the  error 
estimates  of  these  parameters.  If  these  incorporated  values 
are  reasonable,  the  Bayesian  technique  becomes  a powerful 
tool,  particularly  when  the  number  of  observations  is  rela- 
tively small.  The  analyst  should  be  cautioned  however,  against 
using  this  tool  for  a small  number  of  observations  and  then 
attaching  a relatively  large  error  to  the  prior  estimate.  This 
in  effect  attributes  a small  level  of  credence  to  the  prior,  and 
in  the  limit,  will  generate  a solution  for  an  underdetermined 
system.  Obviously,  the  confidence  level  for  the  evolved  para- 
meter estimate  will  be  quite  small  indicating  the  solution  to 
be  useless. 


SUMMARY 


The  different  types  of  optimization  which  have  been  analyzed 
here  can  be  summarized  by  catagorizing  the  estimator  type, 
the  structure  of  the  covariance  matrix  for  the  observations, 
and  the  structure  of  the  covariance  matrix  for  the  prior 
estimates . 

The  estimator  type  can  be  linear  or  non-linear.  When  a 
linear  estimator  is  employed,  no  iterations  are  required. 

For  this  case,  the  general  nonlinear  formulation  can  be 
employed  with  the  initial  estimate  r£  and  the  initial  cal- 
culated response  ug , both  set  equal  to  zero.  From  Eq.  (79) 
the  revised  estimate  is 

r = r + G(u  - Tr  ) (92) 

o o o 

where  G is  defined  by  Eq.  (80).  The  revised  estimate  can  be 
obtained  alternatively  from  Eq.  (77) 

r = (w  + TTwT ) ~ 1 ( w r + TTwu  ).  (93) 

o o 

In  either  case,  the  solution  is  obtained  without  iteration 
by  applying  either  Eq.  (92)  or  (93) . The  revised  covariance 
matrix  is  given  by  Eq.  (83)  . 

The  structure  of  the  covariance  matrix  for  the  observations 
can  be  either  identity,  diagonal,  or  full.  The  case  where 
the  matrix  is  identity  corresponds  to  unweighted  least  squares. 
When  the  matrix  is  non-unit  diagonal,  the  residuals  are 
weighted  with  no  cross-correlation  between  the  observations. 
When  the  covariance  matrix  is  non-diagonal,  cross-correlation 
among  the  observations  has  been  assumed.  As  has  been  noted 
earlier,  the  covariance  matrix  for  the  observations  will  always 
be  symmetric. 
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The  covariance  matrix  for  the  prior  estimate  can  be  either 
infinite,  diagonal  or  full.  The  case  where  the  elements  of 
the  matrix  are  infinite  in  value,  corresponds  to  complete 
lack  of  confidence  in  the  prior  estimate  and  is  equivalent 
to  the  absence  of  a Bayesian  prior.  The  associated  weight 
matrix,  w,  will  be  null,  and  the  estimator  will  attempt  to 
drive  the  calculated  response  toward  the  data.  When  the 
covariance  matrix  is  diagonal,  a prior  estimate  has  been 
specified,  but  no  cross-correlation  is  postulated  among 
its  elements.  Lastly,  when  the  covariance  matrix  is  full, 
cross-correlation  has  been  assumed  between  the  components  of 
the  prior  estimate.  As  with  the  observation  covariance  matrix 
the  covariance  matrix  for  the  prior  estimate  will  always 
be  symmetric. 

The  various  combinations  of  the  two  governing  covariance 
matrices  are  depicted  in  Table  2.  The  four  boxes  on  the  lower 
right-hand  side  give  the  possible  combinations  for  a Bayesian 
prior  estimate  and  are  therefore  separated  from  the  rest 
of  the  table  by  a bold  line.  Note  the  two  empty  boxes  on 
the  left-hand  side.  This  is  because  the  use  of  a diagonal  or 
full  prior  covariance  matrix  along  with  an  identity  observa- 
tion covariance  matrix  would  result  in  a system  of  incompatibl 
units  for  the  objective  function  of  Eq.  (63) , and  is  therefore 
meaningless.  It  should  also  be  noted  that  the  seven  cases 
shown  in  the  table  can  be  applied  with  either  a linear  or 
nonlinear  estimator  giving  a total  of  fourteen  possible  com- 
binations . 

The  revised  estimates  and  covariances  for  the  various  methods 
discussed  here  are  summarized  in  Table  3.  For  the  linear  case 
no  iterations  are  required  and  the  solution  is  immediate.  For 


COVARIANCE  MATRIX  FOR  PRIOR 


COVARIANCE  MATRIX  FOR  OBSERVATIONS 


IDENTITY 

DIAGONAL 

FULL 

INFINITE 

UNWEIGHTED  LEAST 
SQUARES 

WEIGHTED  LEAST 
SQUARES  AND  NO 
PRIOR 

MINIMUM  VARIANCE 

AND  NO  PRIOR 

DIAGONAL 

WEIGHTED  LEAST 
SQUARES  AND  NO 
CROSS-CORRELATION 
ON  PRIOR 

MINIMUM  VARIANCE 
AND  NO  CROSS- 
CORRELATION ON 

PRIOR 

FULL 

WEIGHTED  LEAST 
SQUARES  AND 
CROSS-CORRELATION 
ON  PRIOR 

MINIMUM  VARIANCE 

AND  CROSS - 
CORRELATION  ON 

PRIOR 

Table  2.  Optimization  Schemes  Resulting  From  Various  Combinations 
of  Governing  Covariance  Matrices.  Portion  of  table  with 
bold-face  border  uses  Bayesian  prior  estimate. 
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the  nonlinear  case,  tne  solution  is  obtained  by  iteratively 
replacing  rg  by  r until  the  difference  between  rg  and  r is 
deemed  to  be  sufficiently  small.  For  the  linear  case,  the 
solution  obtained  is  that  for  the  true  extremum  of  the  objective 
function.  For  the  nonlinear  cases,  the  method  may  converge 
to  a local  minimum  or  may  even  diverge.  When  this  occurs, 
special  methods  should  be  employed  to  increase  the  so-called 
radius  of  convergence. 

The  analysis  shown  here  does  not  include  possible  errors  in 
the  measurement  of  the  independent  variables  since  they  are 
usually  relatively  negligible.  However,  there  are  situations 
in  which  both  the  dependent  and  independent  variables  contain 
errors  of  equal  magnitude.  For  these  cases,  a new  estimator 
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NOMENCLATURE 


= matrix  for  linear  set  of  equations 
= determinant  of  Hilbert-like  matrix 
= expectation  operator 
= objective  function 
= linear  least  squares  operator 
= linear  transformation  matrix 
= number  of  independent  variables 
= number  of  observations 
= parameter  vector 

= perturbation  vector  for  nonlinear  estimation 
= covariance  matrix 
= sensitivity  matrix 
= number  of  parameters 
= response  vector 

= residual  or  difference  between  observed  and  calculated 
response 

= free  vector  for  linear  set  of  equations 
= observation  weight  matrix 
= parameter  weight  matrix 
= random  error 


Subscripts 


= latest  iterative  estimate 

= known  prior,  either  experimentally  or  analytically 
= associated  with  parameter 
= associated  with  observation 


Superscripts 


* = revised 
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APPENDIX  B 


VERIFICATION  AND  CHECK-OUT  OF  THE  COMPUTER  PROGRAM,  PEBLS* 

SUMMARY 

When  the  computer  program  "fitted"  Brown  and  Black's  test  data  using 
a near-linear  force-deflection  characteristic,  there  occurred  the  gnawing 
uncertainty,  "Was  the  computer  program  coded  correctly?"  To  verify  that  in- 
deed the  program  was  working  properly,  it  was  decided  to  estimate  the  para- 
meters of  a bi-linear  textbook  example  (for  which  an  analytical  solution  was 
available) . The  estimation  procedure  was  successful  for  this  example,  and 
the  results  are  reported  in  this  appendix. 

B.l  Results  from  Biggs'  Text  (Reference  2) 

There  are  a few  non-linear  problems  of  elasto-plastic  systems  where 
exact  solutions  are  available.  One  such  exact  solution  is  given  in  Biggs' 
text.  Reference  2,  Section  2.7  (p.  69  ff ) . Figure  B-l  describes  the  problem, 
namely  an  elasto-plastic  beam  subjected  to  a Heaviside  step-function  loading. 
Biggs  gives  a graph  of  the  response  (see  Figure  B-l)  and  also  formulas  from 
which  the  response  u(t)  can  be  computed  at  any  time.  (For  further  details 
consult  Reference  2.) 


‘Parameter  Estimation  of  Blast-Loaded  Structures 
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nor**  2 22  Example.  Reeponae  of  elaato-plactir  one-degree  etrueture  to 
suddenly  applied  constant  font. 


Figure  B-l.  Exapmle  From  Biggs  Text 


B . 2 Input  Data  to  PEBLS  Computer  Program 

The  input  data  used  for  this  example  is  shown  below.  Note  that 

the  initial  estimate  of  P. . was  45.3  and  the  initial  estimate  of  d. . was 

lin  lin 

.543,  exactly  the  values  from  Biggs  text  (cf.  Figure  B-l) . These  two  para- 
meters were  input  with  a confidence  of  10%  on  each,  as  shown  in  the  input 
data.  The  "experimental"  data  (which  in  this  case  was  computed,  using  the 
equations  from  Biggs'  text)  is  shown  on  page  B-4.  This  "experimental"  re- 
sponse was  taken  at  seven  time  points  (t^,  t2»  ...t^)  and  was  assumed  to  be 
known  with  an  accuracy  of  about  2%.  There  were  just  two  points  required  to 
define  the  pressure- time  history,  since  it  was  a step-function. 

PARAMETER  ESTIMATION  FOR  BLAST-LOADED  structures 

i-oor  elastic-plastic  test 

I T E RAT I OR  LIMIT  * 12 

BACKSTEP  LIMIT  * 4 

CONVERGENCE  CRITERION  ON  OBJECTIVE  FUNCTION  a 0.00500 
“CONVERGENCE  CRITERION  ON  CHANGE  IN  PARAMETERS  * 0.01000 

MOUSE  STEP  REDUCTION  FACTOR  » 1.00000 

PARAMETER  MOVE  LIMIT  a 0.20000 
THE  TIME  HISTORT  USES  AT  MOST  160  TIME  POINTS 
THE  TIME  INCREMENT  * 0.50000E-03  SECONDS 

THE  MAXIMUM  TIME  THUS  * 0.79500E-01  SECONDS 

THE  SLAB  MASS  PER  UNIT  AREA  * 0.25907E-01 

INITIAL  ESTIMATES  OF  PARAMETERS 


ESTIMATION 

INITIAL 

PERCENT 

PARAMETER 

CODE 

ESTIMATE 

VARIANCE 

CONFIDENCE 

PLIN 

1 

45.300 

20.521 

10.000 

DLIN 

1 

0.54300 

0 . 29485E-02 

10.000 

EPLS 

0 

0.00000 

0.00000 

0.00000 

EMU0 

0 

1.0000 

0.00000 

0.00000 

EMU1 

0 

1.0000 

0.00000 

0.00000 

EMUD 

0 

9.9000 

0.00000 

0.00000 

ALP0 

0 

1.0000 

0.00000 

0.00000 

ALP1 

0 

1.0000 

0.00000 

0.00000 

i ALPO 

0 

9.9000 

0.00000 

0.00000 

BETA 

0 

0.00000 

0.00000 

0.00000 

PSCL 

0 

1.0000 

0.00000 

0.00000 

THE  TOTAL 

NUMBER  OF 

parameters 

BEING  ESTIMATED 

(NP)  ■ 2 
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EXPERIMENTAL  BLAST  AND  RESPONSE  OBSERVATIONS  • DATA  SET  1 
DISPLACENENT  RESPONSE 
“THE  NUMBER'OF  PRESSURE  DATA- POINTS  * 2“ 

THE  NUMBER  OF  RESPONSE  DATA  POINTS  * 7 


EXPERIMENTAL  RESPONSE  DATA 


POINT 

1 

2 

3 

4 

5 _ 

6 
7 


TINE  (SEC) 

0.10000E-01 

0.20000E-01 

0.30000E-01 

0.371C0E-01 

0.45000E-01 

0.55000E-01 

0.66900E-01 


RESPONSE 

0.56530E-01 

0.20836 

0.40781 

0.54300 

0.66363  

0.76352 

0.80551 


VARIANCE 
0.12780E-05 
0.17370E-04 
0.66S20E-04 
0.1 1 790E-03 
0.17620E-03 
0.23320E-03 
0.25950E-03 


CONFIDENCE  (() 
0.19998E-01 
0 . 2 000  3E-01 
0.19999E-01 
0.19997E-01 
0.20002E-01 
0.20001E-01 
0.19999E-01 


BLAST  RECORO 

POINT  TINE  (SEC)  PRESSURE 

1 0 . 00000E  + 00  0 . 30000E+02 

2 0.10000E+01  0, 30000E*02 
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B.3  Initial  PEBLS  Results  for  Biggs  Example 

The  results  of  the  first  iteration  of  the  program  PEBLS  are  shown  in 
Table  B-l.  Note  particularly  that  the  analytic  response  (computed  within 
PEBLS  using  numerical  integration)  did  not  agree  exactly  with  the  observed 
response  (computed  by  hand  using  Biggs  analytical  equations  2.57a,  and  2.61, 
p.  75,  Reference  2).  Also  note  that  the  RMS  error  in  observations  was  com- 
puted as  .889x10  i.e.,  .000889,  which  gives  the  reader  an  idea  of  the 

"goodness  of  fit"  between  the  model  (u-analytic)  and  the  data  (u-observed) . 

After  three  iterations,  the  results  were  as  shown  in  Table  B-2.  Note 
that  now  u-analytic  (computed  by  PEBLS)  is  slightly  closer  to  u-observed 
(the  "experimental"  data)  and  the  RMS  error  in  observations  has  been  halved, 
i.e.,  it  is  now  .399x10  ^ (.000399).  The  final  information  output  by  the 
program  is  shown  in  Table  B-3)  . Note  that  the  parameters  (P^n  and  d^n) 
have  been  estimated  at  45.314  and  .54443,  respectively,  although  Biggs 
exact  values  are  45.3  and  .543.  Thus,  the  estimation  procedure  will  not 
exactly  reproduce  the  parameters  of  our  mathematical  example.  This  (slight) 
discrepancy  may  be  due  to  small  errors  in  numerical  integration  (which, after 
all,  is  not  exact)  or  due  to  the  fact  that  the  "test  data"  were  input  with 
a confidence  of  2%  (and  were  thus  not  known  exactly) . 
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Table  B-l 


computed  and  observed  displacement  response 


TEST  TI Ml 

0.10000E-01 

0.20000E-01 

0, 30000E-01 

0.371001-01 

0.45000E-01 

0.55000E-01 

0.66900E-Q1 


U-ANALYTIC 
0.56367E-01 
0,70780 
0 ,4068  2 
0.542*0 
0.66)10 
0.76250 
0.80)82 


UNOBSERVED 

0.56530E-01 

0.20836 

0.40781 

0.54300 

0.66  36-3 

0.76352 

0.80551 


YaU-UTEST 
•0.16342E-03 
-0.56424E-03 
-0.99456E-03 
-0.19620E-03 
-0.52977E-03 
-0.10177E-02 
-0. 16868E-02 


OBJECTIVE  DUE  TO  OBSERVATIONS 
OBJECTIVE  DUE  TO  PARAMETERS 


0.71 4 206E-01  RMS  ERROR  IN  OBSERVATIONS 
0.000000  TOTAL  OBJECTIVE  FUNCTION 


0.889176E-03 

0.714206E-01 


Table  B-2 


COMPUTED  AND  OBSERVED  DISPLACEMENT  RESPONSE 


POINT 

1 

2 

3 

4 

5 

6 
1 


TEST  TIME 

0. 1 OOOOE-O 1 

0.20000E-01 

0.30000E-01 

0.37100E-01 

0.45000E-01 

0.55000E-01 

0.66900E-01 


U-ANALYTIC 

0 . 56  3 7 0E-01 

0.20785 

0,40706 

0.54331 

0.66396 

0.76376 

0.9054% 


U-OBSEKVED 

O.5653UE-01 

0.20836 

0.40781 

0.54300 

0.66363 

0.76352 

0.80551 


YaU-UTEST 

-0.15988E-03 

-0.51105E-03 

-0.75332E-03 

0.31064E-03 

0.32806E-03 

0.23720E-03 

-0.27962E-04 


OBJECTIVE  DUE  TO  OBSERVATIONS  a 0.452413E-01  RMS  ERROR  IN  OBSERVATIONS 
OBJECTIVE  DUE  TO  PARAMETERS  a 0.706864E-03  TOTAL  OBJECTIVE  FUNCTION 


0.399177E-03 

0.459481E-01 


Table  B-3 


PARAMETER  CONFIDENCE  LEVELS 
INITIAL  LATEST 

PARAMETER  ESTIMATE  ESTIMATE  VARIANCE 

PLIN  45.300  45.314  1.9332 

DLIN  0.54300  0.5444)  0.74290E-03 


PERCENT 

CONFIDENCE 

3.0693 

5.0196 


B.4  PEBLS  Results  for  a Poorer  Initial  Model 

If  the  initial  model  had  been  in  error,  say  one  had  P of  49  and 
d^n  of  .45  (versus  the  exact  values  of  45.3  and  . 54  3 from  Biggs  text),  the 
question  was,  "How  well  will  PEBLS  perform?"  The  input  data  for  this  case 
is  shown  on  page  B-8.  Note  that  now  the  confidence  in  the  parameters  is 
given  as  20%,  whereas  formerly  it  was  10%  (cf.  page  B-3).The  "experimental" 
response  data  were  taken  to  be  the  same  as  previously;  (see  page  B-4) . The 
results  of  the  first  iteration  with  this  input  data  are  shown  in  Table  B-4. 
Note  that  the  maximum  value  of  u-analytic  was  calculated  as  .51185,  whereas 
the  "data"  gave  .80551  for  the  same  time  point,  point  number  7.  Also  note 
that  the  RMS  error  in  observations  went  up  to  .1401,  whereas  in  the  previous 
example  it  was  .000899,  initially. 

After  the  three  iterations,  PE3LS  gave  the  results  shown  in  Table 

B-5.  Note  that  the  RMS  error  has  been  reduced  (from  .1401)  to  .0046,  and 

a significant  improvement  in  "fitting"  the  dynamic  data  had  been  achieved. 

The  final  estimated  values  of  the  parameters  P, . and  d, . are  also  shown 

lin  lin 

in  Table  B-5,  where  the  former  (P^n)  is  given  as  43.8  and  the  latter  (d^n) 
is  .51285.  For  comparison  purposes,  the  "exact"  values  (used  in  Biggs 
example)  were  45.3  and  .543,  respectively.  It  is  noteworthy  that  a rela- 
tively "good"  fit  to  the  data  (RMS  error  of  .0046)  was  obtained,  although 

the  estimated  parameters  p, . and  d. . are  only  within  about  5%  of  their 

lin  lm 

"exact"  values. 


POOR  INITIAL  MODEL 


PARAMETER  ESTIMATION  FOR  BLAST-LOADED  STRUCTURES 


l-Dor 

CLASTIC-PLASTIC  TEST 

ITERATION 

LIMIT  ■ 

12 

BACKSTEP 

LIMIT  * 

4 

CONVERGENCE  CRITERION 

ON  OBJECTIVE  FUNCTION  « 

0.00500 

"CONVERGENCE  CRITERION 

ON  CHANGE  IN  PARAMETERS  ■ 

0.01000 

MOUSE  STEP  REDUCTION 

FACTOR  ■ __  1.00000 

- 

PARAMETER 

MOVE  LIMIT 

* 0.20000 

THE  TIME 

HISTORY  USES 

AT  MOST  160  TIME  POINTS 

THE  TIME 

INCREMENT  * 

0 . 50000E-03  SECONDS 

THE  MAXIMUM  TIME  THUS 

b 0.79500E-01  SECONDS 

THE  SLAB 

MASS  PER  UNIT  AREA  ■ 0.25906E-01 

INITIAL  ESTIMATES  OF  PARAMETERS 

ESTIMATION 

INITIAL 

PERCENT 

parameter 

CODE 

ESTIMATE  VARIANCE 

CONFIDENCE 

PLIN 

1 

49.000  96.040 

20.000 

DLIN 

1 0 

.45000  0 . 8 IOOOE-02 

20.000 

EPLS 

0 0 

.00000  0.00000 

0.00000 

EMUO 

0 

1.0000  0.00000 

0.00000 

EMU1 

0 

1.0000  0.00000 

0.00000 

EMUO 

0 

9.9000  0.00000 

0.00000 

ALPO 

0 

1.0000  0.00000 

0.00000 

ALP  1 

0 

1.0000  0.00000 

0.00000 

ALPD 

0 

9.9000  0.00000 

0.00000 

BETA 

0 0 

.00000  0.00000 

0.00000 

PSCL 

0 

1.0000  0.00000 

0.00000 

THE  TOTAL 

NUMBER  OF  PARAMETERS  BEING  ESTIMATED 

CNP)  B 2 
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Table  B-4 


COMPUTED  AND  OBSERVED  DISPLACEMENT  RESPONSE 
POINT  TEST  TIME  U-ANALITIC  UNOBSERVED  V*U-UTEST 

1 0.10000E-01  0.5S906E-01  0.56530E-01  -0.62388E-03 

2 Q.20000E-01  0.20094  0.20836  -0.74243E-02 

3 0.30000E-01  0.37623  0.40781  -0.31580E-01 

4 0.37100E-01  0.47976  0.S4300  -0.63235E-01 

5 0 , 4 500  0E>  0 1 O.SS280  0.66363  >0.11083 

J 0.5>vvOE>01  0.S7960  0.76362  >0.18392 

7 0.66900E-01  0. 51185  0.80551  -0.29366 

OBJECTIVE  DUE  TO  OBSERVATIONS  • 96.9008  RMS  ERROR  IN  OBSERVATIONS  « 0.140101 

OBJECTIVE  DUE  TO  PARAMETERS  ■ 0.000000  TOTAL  OBJECTIVE  FUNCTION  a 96.9008 


Table  B-5 


computed  and  observed  displacement  response 

POINT  TEST  TIME  U-ANALTTIC  U-UBSEHVED  Y=U-UTEST 

1 0 . 1 OOOOE-O 1 0.S63J1E-U1  0.S6530E-01  >0.19898E>03 

2 0.20000E-01  0.20726  0.20836  -0.11113E-02 

3 0.30000E-01  0,40433  0,40781  -0.34B29E-02 

4 0.37100E-01  0.53761  0.54300  -0.63946E-02 

3 0.45000E-01  0.65640  0.66363  >0.72326£>02 

6 0.55000E-01  0.75901  0.76352  -0.45133E-02 

7 0.66900E-01  0.81159  0.80551  0.60765E-02 

OBJECTIVE  DUE  TO  OBSERVATIONS  a 0.169110  RMS  ERROR  IN  OBSERVATIONS  ■ 0.466160E-02 

OBJECTIVE  OUE  TO  PARAMETERS  ■ 0.767073  TOTAL  OBJECTIVE  FUNCTION  a 0.936183 


PARAMETER  CONFIDENCE  LEVELS 
INITIAL  LATEST  PERCENT 


PARAMETER  ESTIMATE  ESTIMATE  VARIANCE  CONFIDENCE 

PLIN  49.000  43.820  B.1278  5.8182 

DLIN  0.45000  0.51283  0.31881E-02  12.547 


■ 


B.5  Conclusions 

Additonal  work  was  done  to  verify  that  sequential  processing  of  the 
experimental  data  gave  similar  results.  Basically,  the  result  of  this  exer- 
cise was  to  convince  the  authors  that  the  program  PEBLS  was  functioning 
correctly  and  numerically  integrating  properly.  Since  exact  analytical  exam- 
ples are  not  available  with  varying  p,  a,  etc.,  all  the  loops  and  logic  of 
PEBLS  have  not  been  verified  by  hand  calculations.  Nevertheless,  the  authors 
can  state  with  some  confidence  that  the  program  is  completely  functional, 
and  it  appears  to  be  properly  coded. 
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